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ABSTRACT 



Advances in digital computer technology have led to employment of 
the computer as a control device. Although many theories for utilization 
of the digital computer as the compensator in a position feedback control 
system have been proposed, there is little documentation of digital 
control performance in a hardware environment. This thesis is an 
investigation of the performance of a real system under sampled-data 
control. The variable gain amplifier sampled-data control theory is 
first simulated and then tested on a real system. The hardware and 
software developments necessary for implementing this control theory 
are discussed in detail. Additionally, this method is incorporated in 
a time-sharing computer program for controlling many systems simultane- 
ously from a single computer. Results of the tests are presented and 
evaluated . 

The authors wish to express their appreciation to Professor Robert 
D. Strum of the U # S. Naval Postgraduate School for his guidance and 
assistance during this investigation. 
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CHAPTER I 



1.1 Introduction 

During the past 15 years, advances in digital computer technology 
have developed the digital computer into an invaluable tool for all 
branches of the scientific world. In almost every scientific discipline, 
the capabilities of the digital computer for vast computational tasks, 
data processing, and simulation have been profitably employed. In 
addition to these uses, some engineers in the automatic control field 
have considered the digital computer, operating in a real time mode, as 
a tool for improving control techniques and capabilities. It is in this 
area of computer application that the authors are extremely interested. 

There are many theoretical expositions on sampled data control 
techniques, and some include results from the digital control of 
simulated plants. However, with the exception of the technical reports 
of the Navy Electronics Laboratory [1, 2,] the authors were unable to 
find any information on sampled-data control of an actual plant. It is 
the opinion of the authors that unless the theory is tested in a hardware 
environment, the results are academic because of the difficulty in 
accurately simulating the inherent non-linearities in a real system. It 
was decided at the outset, therefore, that a prime goal of the thesis 
was to test the theories of sampled-data control on a real plant. 

A second area in which data was lacking involved one of the primary 
advantages claimed for the sampled-data control techniques: the cap- 

ability to control many systems simultaneously from the same computer. 
Although implementing such a multi -control system is primarily a computer 
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programming task, certain basic techniques must be developed to insure 
flexibility and adequate performance. It was felt that the lack of 
documentation in this area was a serious shortcoming. Thus, a second 
goal of this thesis was to develop and implement techniques for multi- 
system sampled-data control. 

Having established the basic goals of the thesis, the authors made 
a decision which essentially limited the scope of the thesis. Throughout 
the literature, there are many proposed theories for sampled-data control 
systems. It was decided to select only one of these theories as the 
basis for the laboratory experimentation. Rather than testing many 
theories in a less rigorous fashion, the authors preferred to make 
exhaustive hardware tests on the basis of one theoretical solution to 
the control problem. In general, the difficulties encountered in one 
theoretical approach would be comparable to those of another, and the 
techniques developed may be extended to other theories. 

After examining many of the proposed theories, the authors, in 
conjunction with their advisor. Professor Robert Strum, selected the 
theory proposed by Benjamin C. Kuo. The basis of this theory is the 
concept of a variable gain amplifier in the control loop. The theory 
provides for minimum time, ripple free response to deterministic inputs. 
The prime reasons for selection of this theory as a basis for the 
laboratory experimentation are the theoretical advantages it offers, 
and the relatively simple tasks the digital computer is required to 
perform. 

In considering the laboratory portion of the thesis, the authors 
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had one more basic decision to make. In selecting a digital computer 
for the experiments, the facilities at the Postgraduate School offered 
a choice between a general purpose computer, the Control Data Corporation 
1604, and a small data processing computer with an auxiliary arithmetic 
capability, the CDC-160 in conjunction with the CDC-168 arithmetic unit. 
Due to the fact that any practical control method must not only perform 
well, but also must be economically feasible, the smaller, less expensive 
160 computer was selected. It is conceded that the larger computer could 
more efficiently perform the required controlling tasks, but the authors 
preferred the more realistic practical approach. This led to a secondary 
goal of the thesis, to be able to comment on the general feasibility of 
sampled-data control systems on the basis of the results obtained. 

1.2 Digital Controller Design 

Before proceeding with work on the real plant, the authors decided 
to select a technique for the design of digital controllers for sampled- 
data control systems. A design theory that would give accurate results 
and would be readily adaptable to the existing hardware was desired. It 
was determined that the digital controller design described in Section 
9-9 of Analysis and Synthesis of Sampled-Data Control Systems by 
Benjamin C. Kuo would be well suited to this problem. [^3^] Kuo’s design 
theory makes use of state space and state transition techniques. It 
also enables derivation of a D(Z) that will yield a minimum time, ripple 
free response to the designated deterministic input (i.e., ’’deadbeat 
response”). This approach to digital controller design describes the 
digital controller as a variable gain amplifier which outputs certain 
values each sampling period in order to control the continuous system 
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in the desired manner. It is easily adaptable to implementation on a 
digital computer, and the use of state variables makes the manipulation 
of the mathematics involved relatively easy. 

To explain this "variable gain amplifier" design theory it is best 
to look at a sample problem. Consider the sampled-data control system 
of Figure 1. The variables and x^ are the state variables for the 
controlled process where x^ = x^ ; h(t) is the output of the zero-order 
hold and e(t) is the actuating signal. Note: e(t) = r (t ) -c (t) ;x^=c (t ) 




D(Z) 



Figure 1 

To determine a D(Z) for "deadbeat response" to a given input: 

(1) Draw the state transition flow graph of the system with the 
digital controller represented by a variable gain amplifier 
K(kT) . K(kT) = h (kT) / e(kT). 

(2) For a "deadbeat response" the system error must be zero for 
t > nT where n is the smallest possible integer. 



i.e. , 


x 1 (nT) = r (nT) 






(nT) = x^(nT) = . . . .= x p (nT) = 0. 


unit step 




x 2 (nT) = 1; x 3 (nT)=x 4 (nT)=. . ,=x p (nT) = 0. 


unit ramp 



From this information the variable gains (K^ 1 s) and the minimum 
n can be determined. 
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(3) Compute the h (0 + ) , h(T + ), . ...h(nT + ) using the 1 s 
determined and recalling that e (kT + )=r (kT) -x^ (kT) . 

(4) Calculate D(Z) using the relationship D(Z)=H(Z) /E (Z) . 
Returning to the example, let us draw the state transition flow graph 
of the system. 




Figure 2 

Using Figure 2 and applying flow graph techniques and the inverse 
Laplace transformation we arrive at the following: 

x^ j^(k+l)T J = jl-K^T-l+e'Lj x ;L (kT)+(l-e" T )x 2 (kT)+K k (T-l+e" T )r (kT) 
x 2 [(k+l)lj = -K k (l-e‘ T )x 1 (kT)+e' T x 2 (kT)+K k (l-e' T )r(kT) 



Now let us assume a unit step input and a sampling period of one tenth 
of a second. ( r(kT) = 1.0, T=0.1) 

X x £ (k-KL)T J = (1 - 0.005K k )x 1 (kT) + 0.095 x 2 <kT) + 0.005^ 
x 2 I (k+l)T 1 = -0.095K k x 1 (kT)+ 0.905x 2 (kT) + 0.095*^ 



For k=0 and (0) = x 2 (0) = 0. 



5 



x - (T) = 0.005K 
1 o 

x 2 (T) = 0.095K q 



From the constraints placed on the problem (i.e., x^ (nT) = r(nT) and 
X 2 (nT)= 0 .) it may be seen that the solution can not be obtained from 
the above. In practice this means that the output of the controlled 
process cannot be made equal to a step input in one sampling period. 
We must therefore proceed to the second sampling period and let k=l. 
Doing this and setting x^(2T)= 1.0 and x^(2T) - 0. yields: 

X, (2T) = 0.014K - 0 . 000025K K. + 0.005K, = 1.0 

1 o o 1 1 

x„(2T) = 0.086K - 0 . 000475K K. + 0.095K.. = 0.0 

2 o o 1 1 

from which 



K = 105.6 ; K, = -203.0 
o 1 

therefore 



e(0 + ) = r(0)-x 1 (0)=l 



e(T + )=r(T)-x 1 (T) = 0.472 



h(0 + ) = 105.6 



h(T + ) = K 1 e(T + ) = -95.8 



and 



D (Z) = 



h(0 + ) + h(T + )z" 1 

e(0 + ) + e(T + )Z _1 



105.6 - 95.8 Z 

1 + 0.472 Z" 



The above example demonstrated the technique of designing a digital 
controller for a sampled-data control system using the variable gain 
amplifier concept. It should be observed that in implementing this 
controller using a digital computer, the D(Z) does not have to be actually 
calculated. In fact, only the K, 1 s need be entered in the computer. 
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This proves to be quite advantageous since for a given sampling period 
the K^'s are only a function of the continuous system gain and the G(S). 
Thus for a digital computer controlling a number of plants the K^'s 
could be pre-computed and stored in the computer ready for use. These 
variable gains when multiplied by the sensed error and output at the 
sampling rate will provide the desired system response to the designated 
input. In fact, the digital controller designed above will yield the 
desired "deadbeat response" to any step input. It must also be pointed 
out that the design given is only for a step input and will not give 
optimum response for a ramp or parabolic input. 

The design problem as stated and solved above was for the control 
of a type one system in response to a step input. In studying this 
design technique a digital controller was also designed for the same 
type one system (G(S)= ^ ^ f° r "deadbeat response" to a unit ramp 

input. With respect to a ramp input, optimum response was obtained in 
three sampling periods. It was also noted that once zero error had 
been achieved a constant output of 1.0 was needed from the controller. 
This was required because a type one system has an inherent steady state 
error in response to a ramp input. In this context it is important to 
realize that with the use of a digital controller a continuous system 
with an inherent steady state error can be made to respond with "deadbeat 
response" and no steady state error to any particular deterministic 
input. (See Appendix I for details.) 

In Appendix I a digital controller is also designed for a type 

1 



zero system (G(S) = 



(s + 1) (s + 2) 



) which has an inherent steady state 
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error to both step and ramp inputs. With the designed controller, these 
errors were reduced to zero. 

With regard to type two systems some difficulty was encountered. 

For a system of the type G(S) = — a digital controller was readily 

S 1 

designed; however, for a system of the type G(S) = — r the mathe- 

s (s+1) 

matics involved in evaluating the K^'s proved to be cumbersome. An 
iterative solution to the non-linear simultaneous equations was used 
in evaluating the K^'s. This proved to be feasible. Further investiga- 
tion into the problem was not carried out since the authors wished to 
direct their efforts to the control of the real system. 

1.3 Simulation and Results 

After completing the design of a digital controller for a "dead- 
beat response" to a unit step input, it was decided to simulate the 
system response. This was done to provide a check on the computations 
and to verify the design theory itself. The simulation programs used 
were written in Fortran 60 and employ a library routine known as INTEG 1. 
INTEG 1, written by Dr. J. R. Ward, provides a fourth order Runge-Kutta 
solution to ordinary differential equations and was well suited to the 
simulation problem encountered. (Appendix II contains a complete 
simulation program.) 

The graphical results to the simulation are on the following pages 
and do verify the design theory as well as the computations involved. 
(Appendix I contains graphical results for a unit ramp input as well 
as results for a type zero system in response to both unit step and ramp 
inputs . ) 
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It should be mentioned that a certain degree of accuracy is needed 
in the determination of the h f s used. In some simulations, the rounding 
of the h value to two decimal places resulted in a response that was 
slightly less than optimum. 



9 



009 092. 991 096 00b 010 012 OH 




000. Ml Ml 003 004 0051 



K-dCft LE ' JJ10E-0J UNJT3< / INCH. 

'f-XALE - O,00E-0J UNJiTXJNCH. 

BROWNE DIGIT 1 

RUN 1 OUTPUT US. TIME 

Figure 3 
i n 



002 001 009 OO'h 010, 012 0 lit 




KSCA LE - J.ME-flJ UNJT3/JHCH. 
'fOMLE ’ Q.ME+M UNJTV1NCH. 



BRQkJNE DIGIT 1 

RUN 1 VELOCITY US.TIME 



Figure 4 
11 




KnSCflLE -v d.ME-«J LJNJTS/JNCH. 
'T-3CALE - UNIT 3 'INCH. 



BRCVJNE DIGIT I 

Ru ^ 1 EDOT ^ ERRDR 



Figure 5 
12 



CHAPTER II 



2.0 Investigation 

The investigation conducted during the course of this thesis was 
developmental in nature. Using the variable gain amplifier control 
theory discussed in Section 1.2, techniques for accomplishing the 
sampled-data control of a real plant were developed. The investigation 
also included the development of a computer time sharing system through 
which many systems can be simultaneously controlled. 

2.1 Theoretical Design 




D(Z) 

Figure 6 



Figure 6 shows the general block diagram of the sampled-data 
control system that was used. The system breaks down into two main 
components which are the continuous plant and the D(Z) , sampler, and 
zero-order hold. The continuous plant is described in Appendix III. 
The CDC-160, CDC-168, and A/D, D/A converters function as the sampler, 
D(Z), and zero-order hold. This operation will be treated in detail 
in Section 2.2.1. 

The digital controller, D(Z), was designed to give a "deadbeat 
response" to a step input using the theory described in Section 1.2. 
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It was decided to solve for the K^'s in general terms initially and 

then substitute the known G(S) determined in Appendix III. 

A 

Figure 7 is shown for a generalized G(S) = s ( s+ ^ ) 




From the flow graph we have 
KA 



, x r i ka i 

‘l (s)= - - 2 ~ ,f l 

L s (s+c*)J 



(t o )+ s(s+o<) ~2 k ^o j 



Xr>(t) + 



KA 



s ( s+o< ) 



, X -KA ^ x . 1 x , KA 

X 0 (s)= — — 7- X- (t ) + — X 0 (t ) + — — r- 

2 s(sH-ck) 1 O S+cX 2 O S ( S+ o< ) 



r(t o ) 



r < t o ) 



and after taking the inverse Laplace transformation we arrive at 

- ocT 



(kT) + (l-e' wT )x 2 (kT)+ 



|(k+l)T]= |l - («T-l+e"°' T )jx 1 

jo^H-l ) tJ = (l-e‘ ° <T )x 1 (kT)+e" WT x 2 (kT)+ -|p (1-e * T ) r(kT) 



cxT-l+e 



- <*T 



r (kT) 



For x^(0) = x^(0) = 0.0 

and for a step input of r(0) = r(T) = (2T) - r 

we have 



K = 
o 



h = 
o 



ex 



- o< T 
AT(l-e ) 



K, = 



.2 -<XT 
CX e 



A [e _ wT (l+«T)-l| 



cX 



AT ( 1-e” * T ) 



Ar 



h l = 



-cXe 



- <*T 



_ at T 

AT ( 1 - e ) 



A r 
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where (2T) = r: x^(2T) = 0.0 

It is seen that a "delta" term is present in the solution for the 
K^'s and the h's. This "delta" term arises from the fact that the 
error signal was attenuated by a 0.0121 factor before being sent to 
the A/D converter. This will be discussed in detail in Section 2.2.2, 



Now for the G(S) = 



40 



s (s-f3) 



; T = 0.3; where r = 5.0 we arrive at: 



K = 0.421 

o 

h = 0.0252 

o 



K x = -0.398 
h = -0.0102 



The h's computed were used in an earlier trial program. The computed 
K^'s were used in the final control program and as such were entered 
directly in the computer. The CDC-160 together with the CDC-168 were 
programmed to serve as the digital controller using these K^ 1 s . The 
CDC-160 also operated as the sampler since the sampling period was also 
programmed into the computer. 

The real system response was also simulated using the same 
simulation scheme as Section 1.3. The simulation results are in 
Section 2.3. 

After controlling the real plant in a "deadbeat response" manner, 

it was decided to work on a time sharing routine for optimum control 

of two plants. The second plant was an analog simulation set up on 

1 



the EAI TR-20 analog computer with a G(S) equal to 



The CDC- 



s (s+1) 

160 was then used to provide optimum control to both systems. The 
design of the digital controller for the second plant followed the same 
theory used before. 
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2.2 Development of Sampled-Data Control Technique 



The step from the theoretical control solution to a working sampled- 
data controller involves the development of methods through which the 
requirements of the theoretical solution may be implemented in a real 
controller. This is a two- fold development involving the software of 
computer programming and the associated hardware requirements. Develop- 
ment in each area will be discussed in detail. 

2.2.1 Hardware Development 

Having developed, theoretically, the performance of the D(Z) 
portion of the control system block diagram, it was necessary to develop 
a hardware equivalent to the previously described mathematical model 
of D(Z). To accomplish this hardware development, the following equip- 
ment was used: 

(1) A two-channel digital to analog converter. (D/A) 

(2) A four-channel analog to digital converter. (A/D) 

(3) A Control Data Corporation 160 computer, in conjunction 
with a Control Data Corporation 168 arithmetic unit. 

(4) Operational amplifiers used for summing of signals coming 
in and out of the converter units. 

The equipment used was chosen, not because it was thought to be most 
efficient for the desired tasks, but because it was readily available 
in the Digital Control Laboratory of the U„ S. Naval Postgraduate School. 
Comments on the relative merits of the equipment in performing these 
tasks will be included where appropriate. 

Development of the hardware equivalent of the mathematical model 
of D(Z) will now be discussed. It was decided immediately that the need 
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for analog hold circuits could be eliminated by replacing them with 
more accurate computer programmed holds. This is possible because of 
the operating characteristics of the converter units in conjunction 
with the 160 computer. First, the 160 computer can sample a converted 
digital input only when specifically commanded, and then it samples 
only one value. Similarly, the 160 computer can output a digital value 
for conversion only upon command, and furthermore, that value remains 
constant in the output register until it is changed by inserting a new 
value. Thus, by programming the computer to sample converted analog 
inputs and provide digital outputs for conversion at specified time 
intervals, hold circuits become unnecessary. In fact, the holds be- 
come an inherent part of the sampling rate of the computer and are 
varied automatically as the sampling rate changes. 

Examining Figure 6, one sees that the input to D(Z) is the error 
signal, e(t), which is the sum of the reference signal, r(t), and the 
negative control signal, -c(t). There are two methods by which e(t) 
may be transmitted to the computer. The first method is to sample r(t) 
on one channel of the A/D and to sample -c(t) on another channel. The 
two inputs could then be summed in the computer to obtain the error 
signal. This method was discarded for two reasons. First, since it 
is impossible to sample r(t) and -c(t) simultaneously, the error signal 
computed would have some inherent error. Second, the A/D conversion 
time is of the order of 120 microseconds, and it was felt that this 
additional time delay might become prohibitive at sampling rates approach- 
ing 0.1 seconds. 
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The second method of supplying the error signal to the computer, 
and the one which was adopted, involved summing r(t) and -c(t) in an 
operational amplifier prior to transmission to the A/D. This method 
has the advantage of being fast, accurate, and readily adaptable to a 
restriction imposed by the physical characteristics of the converter 
unit. This restriction is that the converter units accept only negative 
voltages in the range of zero to minus ten volts. Thus, to handle both 
positive and negative values of e(t), the sum of r(t) and -c(t) is 
added by means of an operational amplifier to minus five volts. A 
zero error signal would then be sampled by the computer as an input of 
minus five volts, and an error signal of three volts would be sampled 
by the computer as an input of minus two volts. Similarly, the 
converted D/A value which is transmitted to the amplifiers of the 
system must be added to five volts prior to insertion in the system. 



A block diagram of the hardware system used to provide the 
equivalent to the mathematical model of D(Z) is shown below: 




Figure 8 



With the above hardware system established, the rest of the development 
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of a unit equivalent to the mathematical model of D(Z) must be 
accomplished by programming the computer. 

2,2.2 Software Development 

Programming the computer to accomplish the aforementioned task 
was done in a progressive fashion, starting with programs to accomplish 
the simplest control functions and expanding the basic programs to 
provide for accomplishing more complex tasks. This was done for two 
reasons. The primary reason was that at the outset the authors wanted 
to isolate quickly any faults in the digital control system, and to 
analyze carefully system response at each stage of control function 
complexity. A secondary reason for adopting this progressive program- 
ming technique was the fact that the 160 computer is rather limited in 
its capability^ . Due to this limitation, the authors desired to bring 
the control capabilities of the equipment used to a maximum before 
computer saturation problems were encountered. Following is a discussion 
of these progressive programs. 

The first program written was intended primarily to check the 
operation of the hardware equivalent of the D(Z) under minimum capability 
requirements, and also to check the plant response with the D(Z) in the 

^“The 160 has limited storage capacity (4K) , and a relatively slow 
memory cycle (6.4 microseconds). It can only multiply or divide in 
conjunction with the 168 and the 900 microsecond time for such operations 
is quite slow. 
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loop. The essence of Program 1 is to enable the plant to run as a 
continuous system with the D(Z) block in the loop. Mathematically, 
such a system may be described as D(Z) = 1., and the sampling period, 
T, approaching zero (approximately 200 microseconds). 

The programming aspects of Program 1 were quite simple. The 
program was an iterative loop in which the current error signal was 
sampled from the A/D converter, read into the computer, stored 
temporarily, and output to the D/A converter for return to the system. 
A flow chart for Program 1 is shown below. A complete text of Program 
1 is presented in Appendix VI. 



Program 1 
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The results of the test using Program 1 were quite satisfactory. 
System responses were obtained from step, ramp, and sinusoidal inputs, 
and these responses matched those obtained from the system when operated 
as a purely continuous system with unity feedback. This, of course, was 
predictable with D(Z) = 1 . , and T approaching zero. However, the results 
do validate the selection of hardware used to implement the realization 
of the mathematical D(Z). 

The results of Program 1 opened another avenue of approach to the 
problem of implementing the variable gain samplifier method of digital 
control. It was noted, when working with the system in the Digital 
Control Laboratory, that the system was subject to many random inputs 
of small magnitude from various sources of noise. Program 1, when 
tested, maintained stability, and accurate and fast response in the 
presence of these noise sources. It was further noted that the advantages 
of the variable gain amplifier method in regard to minimum time response 
are lessened when the inputs are small and random. 

In the light of the foregoing facts, a temporary control philosophy 
was developed. This philosophy states that for step inputs which exceed 
in magnitude a threshold, a standard variable gain amplifier control 
solution would be employed. For small inputs, such as noise, the system 
would be controlled as it was with Program 1. The threshold would be 
set at a level which would prevent the implementation of a solution of 
the type described in Section 1.2 for the random noise inputs, which 
were restricted to a relatively low magnitude. In other words, the 
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system would run in an essentially continuous mode D(Z)=1., until 
receipt of a valid step input, at which time it would be controlled 
by a variable gain amplifier digital solution. 

To implement this type of control philosophy, an intermediate 
program, Program 2, was written to modify Program 1 to provide for 
detection of a definite step input and an exit to a solution. This 
detection is accomplished by comparing the new value of the sampled 
error signal with the immediately preceding value of the sampled error 
signal. If the difference between these values is greater than the 
established threshold, a step input is detected. To assist the reader 
in understanding the programming procedure for determining the sign of 
a detected step input, the following table of octal number values for 
analog voltages into the converter is presented: 



ANALOG OCTAL 

VOLTAGE NUMBER 



0 


4000 


-1 


4632 


-2 


5463 


-3 


6314 


-4 


7144 


-5 


0000 


-6 


0631 


-7 


1463 


-8 


2314 


-9 


3144 


-10 


3777 



A flow chart for Program 2 is presented on the following page. A 
complete text of Program 2 is presented in Appendix VI. 

The results of Program 2 were excellent. It proved to be a fast 
and accurate method of sensing step inputs having a magnitude greater 
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Program 2 




Figure 10 
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than the threshold level. In one sense Program 2 was too sensitive* 

When certain knife switches were used to provide zero to five volt step 
inputs, negative transients appeared at the instant the switch was thrown. 
Program 2 was so sensitive that it sensed the transient as a negative 
step input. This imposed the requirement that step inputs be obtained 
from an electronic switching device or from a high quality knife switch 
which eliminates the aforementioned transient. 

The next step in the development of the temporary control philosophy 
involved utilizing the computer in a real time mode to implement a 
variable gain amplifier solution. The iterative loop of Program 2, 
which maintains continuous control until a step input is sensed, employs 
the computer in a free-running mode. The variable gain amplifier solution, 
however, requires that the computer insert control voltages of definite 
real time length into the system. To do this, the computer must output 
a control voltage to the system, hold the voltage for a specified length 
of time, and then output another control voltage. Program 3 was written 
to develop techniques for accomplishing this task. 

Program 3 was a simple program which called for the computer to 
output a given voltage, delay through a timing chain, and output another 
voltage. Theoretically, the length of the delay could be calculated 
from the execution time for the instructions which constitute the timing 
chain. However, exact timing data for D/A conversion were not available. 
Therefore, approximate timing chain delay was calculated and adjusted 
by experimental results for exactness. The timing chain was formed by 
constructing an iterative loop in which an index was increased by one 
on each pass and compared to a preset total. When the index equaled the 
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preset total, the time delay was complete. A flow chart for Program 3 
is shown on the following page. A complete text of Program 3 is 
presented in Appendix VI. 

In testing Program 3, it was decided to design a basic delay block 
of 0.1 second. With an accurate delay block established, delays which 
were integer multiples of 0.1 second could be achieved by iterating 
through the basic delay the required number of times. Delays of less 
than 0.1 second could be achieved by setting the preset total at a proper 
fraction of that required for a 0.1 second delay. 

Satisfactory results were obtained from Program 3 in that a delay 

of exactly 0.1 second was produced between the output of the first and 

second voltages. For this delay, the preset total was set at 3403 o . 

o 

Results were checked on a high speed Mark II Brush Recorder and an 
oscilloscope . 

Having developed Programs 2 and 3, the basic tools for implementing 
a solution within the constraints of the temporary control philosophy 
were ready. By using Program 2 to maintain continuous type control in 
the absence of an input signal and to sense a step input, and by employing 
the features of Program 3 to output the required controlling voltages, 
the desired solution to a step input may be attained. To test the basic 
validity of their approach, the authors wrote Program 4 as a first trial 
of the variable gain amplifier control method. 

Program 4 was an unsophisticated approach which started with the 
system being controlled in the undisturbed state by Program 2 which will 
hereafter be referred to as the Sense Loop. When an input is sensed, h^ 
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Program 3 




Figure 11 

is transmitted immediately to the system. is held for one sampling 

period, and then h^ is transmitted and held for one sampling period. At 
the completion of the second sampling period, the system returns to 
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continuous operation as the computer returns to the Sense Loop. To 
insure a return to the continuous mode of operation at the completion 
of the second sampling period, it was necessary to "bootstrap" the Sense 
Loop by sensing the current error and inserting it in the Sense Loop as 
the M old value". Were this procedure not followed, an imperfect solution, 
one which did not have zero error at the end of two periods, would cause 
the Sense Loop to sense another step input and thus cause the system 
to go into unstable, erratic operation. 

The unsophisticated aspect of Program 4 stems from the fact that 
the h’s were not computed by multiplying e(t) by K during the solution, 
but were pre-calculated and inserted in memory as constants to be trans- 
mitted to the system at the proper time. Kq and have been determined 
in Section 2.1. h^ was simply times e^(t), and h^ was times e^(t) 
(predicted). Due to the fact that r(t) and -c(t) were transmitted to 
the summing operational amplifier after passing through dropping resistors 
in the front end of the system, the magnitude of the error signal was 
considerably decreased. For a five volt step input, e^t) sensed by 
the computer was 0.0600 volts and predicted e^(t) was 0.0258 volts. On 
the basis of these values, h^ = 0.0252 volts and h^ = -0.0102 volts. A 
0.3 second sampling interval was used. 

A flow chart for Program 4 is shown on the following page. A 
complete text of Program 4 is presented in Appendix VI. 

Complete results of the responses obtained when using Program 4 
will not be presented because these tests were in the nature of feasibil- 
ity checks. The results may be summarized by stating that they provided 
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Program 4 




Figure 12 
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definite proof of the validity of the variable gain amplifier approach, 
and led the authors to adopt certain techniques which were used in 
later programs to improve performance. Some of these lessons learned 
from these tests will now be discussed. 

Accuracy problems were caused by the small magnitude of the error 
signal received at the D(Z) portion of the system for a five volt 
step input. Due to this fact, the magnitude of h's to be transmitted 
to the first amplifier was so small that the resolution of the D/A 
converter prevented sufficient accuracy for these h's. The resolution 
for the converter is 0.0024 volts per octal number. With h's of the 
magnitude of 0.02 volts, accuracy was limited to the first digit of 
the h. This condition was barely satisfactory due to the stringent 
demands of the variable gain amplifier method for accuracy in the h's. 
This demand was demonstrated when solutions were simulated on a 
general purpose digital computer. This demand was discussed more 
thoroughly in Section 1.3. 

To correct the problems associated with converter resolution, 
the obvious solution was to increase the voltage levels in and out 
of the converter. The simplest method of accomplishing this would 
have been to eliminate the dropping resistors before the summing 
amplifiers. However, when this was tried, system noise increased to 
a level which was not tolerable. The method of achieving increased 
accuracy, which produced the best results, was to multiply the error 
signal by ten in the summing operational amplifiers prior to A/D 
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conversion, and to divide the D/A converted voltage by ten in the 
summing operational amplifiers before transmittal to the system. This 
enabled the use of h's of ten times computer value for conversion 
which alleviated the resolution problem. This scheme also maintains 
operation in the continuous mode of the Sense Loop to continue as 
before. For all further tests, the hardware shown in Figure 8 was 
modified to provide for the proper multiplication and division in 
the summing operational amplifiers. 

A second lesson learned in this test program related to the 
"bootstrap" procedure discussed previously. This feature was 
originally inserted in the program to prevent system runaway for 
incorrect solutions resulting from tests using improper h's. It was 
found, however, that this "bootstrapping" return to continuous operation 
at the completion of the second period was an integral part of a 
satisfactory solution. This is due to the fact that a real system 
has inherent non-linearities such as coulomb friction and backlash. 

In view pf these non-linearities, a perfect theoretical solution can 
not be obtained. (It is possible to approach more closely the perfect 
solution by slight modification of the size of the h's.) Thus, a 
rapid, smooth return to the continuous mode after the second sampling 
period is essential to a satisfactory solution. It is noted that 
by properly adjusting the size of the h's, the error at the end of 
the second sampling period is small (less than 10% of the step size), 
and the use of the continuous mode reduces error to zero rapidly. 
Specific examples of the foregoing may be noted in the presentation 
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of results in Section 2.3. At that time, the reader may note that 
the necessity to go to the continuous mode to bring error to zero 
does not significantly degrade system performance in comparison to a 
theoretical solution . 

With the information derived from Programs 1 through 4, the 
authors developed Program 5 which was intended to culminate efforts 
in the development of the temporary control philosophy. Program 5 was 
a generalized version of Program 4, the primary difference being that 
in the newest program the h's were computed during the solution 
instead of using precomputed values. The solution technique may be 
summarized by stating that each h is computed by sampling current 
error and multiplying by the appropriate gain constant (K^) by using 
the CDC-168 arithmetic unit in the multiply integer mode in conjunction 
with the CDC-160 computer. 

Programming arithmetic operations for these two units is quite 
straightforward when Subroutine Arith, written by Professor M # L # 
Cotton, is used. However, the characteristics of the CDC-168 required 
manipulation of the sensed error signals. Specifically, for an 
arithmetic operation, the 168 requires two 22-bit operands, and 
supplies a 22-bit solution. Each operand is composed of two 12-bit 
160 words, with the least significant half of the number in an even 
numbered cell and the most significant half of the number in its odd 
numbered mate. The first bit in each cell pair is a sign bit. It is 
this sign bit in the first bit of each of the two words composing the 
operands that require the manipulation. For the K multipliers which 
remain constant, proper values may be inserted into both storage cells 
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manually before the program runs. However, for the sensed error 
signal, negative numbers received from the converter create a problem. 
For the 168 to function properly, the sign bits for both portions of 
the operands must agree. To handle this problem, the following proce- 
dure was employed. To derive the proper h, the product of and e(t) 
is required. After e(t) is sensed, it is placed in cell 0024 as the 
least significant half of an operand. The most significant half, 
cell 0025, is preset to zero. For negative step inputs, the error 
signals read in from the A/D converter are positive numbers, and the 
multiplication process may proceed immediately. For positive step 
inputs, however, the error signals sensed from the converter are all 
negative numbers. To maintain sign bit consistency with cell 0025, 
this negative number is complemented and placed in cell 0024. The 
multiplication operation then proceeds. The product, of course, is 
of the wrong sign so it must be complemented prior to being transmitted 
to the system. This procedure will be defined completely in the flow 
chart for Program 5. 

The other major innovations of Program 5 involved attempts to 
obtain greater accuracy in the solution ^nd maintain high resolution 
in the converter. As previously discussed, the multiplication and 
division by ten in the operational amplifiers were included. To 
achieve greater accuracy in computing the h's, 100 times was set 
in the computer as the multiplying constant. The products were then 
divided by 100 prior to being output to the D/A. 
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As before, a sampling interval of T = 0.3 seconds was used. 

The gain constants for this system were computed to be: Kq= 0.421 and 

= -0.398. A flow chart for Program 5 is shown below. A complete 
text of Program 5 is presented in Appendix VI. 



Program 5 
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As noted before, Program 5 marked the culmination of efforts 



in conjunction with the temporary control philosophy. Although the 
results of Program 5 were good, and gave valuable insight into digital 
control technique, Program 5 was merely a vehicle through which the 
validity of variable gain amplifier theory could be tested. The 
temporary control philosophy, although perfectly satisfactory for 
control of a single system, is unsatisfactory as a general philosophy 
for digital control because it excludes utilization of one of the 
primary advantages claimed for digital controllers: simultaneous 

control of many systems by the same computer. The authors believe 
that the realistic application of digital control methods is closely 
tied to the capability for simultaneous control of many systems. 

Toward this end, a final control philosophy was developed. 

The goal of this final control philosophy was to achieve this 
simultaneous control. This goal implies that the computer be "time- 
shared 11 among the systems to be controlled. Therefore, the significant 
difference between the temporary and final control philosophies is 
one which enables the computer to time-share its control function. 

Before discussing the final control philosophy in detail, 
specific performance goals for this philosophy will be discussed. 

The authors wanted to develop techniques for controlling two or more 
systems simultaneously by employing one digital computer on a time- 
share basis. The control was to be achieved by the variable gain 
amplifier method. To achieve generality, it was desired to include 
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the capability for the systems to be controlled at different 
sampling rates. Furthermore, it was considered necessary to maintain 
the capability to handle simultaneous step inputs to each system 
without performance degradation. 

It may be recalled that the temporary control philosophy called 
for a return to a f ree-running Continuous 11 mode at the completion of 
the solution. The major difference in the final control philosophy 
is that at the completion of the second sampling period, the computer 
returns to a mode in which D(Z) is still 1.0, but instead of a sampling 
interval approaching zero, a sampling period of definite finite 
length is employed. Although this innovation does not change the 
response characteristics of the control solution, it marks the change 
from a partially discrete, partially continuous system, to one which 
is completely discrete. It is this change which enables the develop- 
ment of the required time-sharing computer operation for multiple 
control . 

Having discussed the major difference in the control philosophies, 
the final control philosophy will be covered in detail. As before, 
with no signal input to either system, the computer resides in a 
Sense Loop. In the Sense Loop each system under the control of the 
computer is sampled and tested for a step input. If a step input 
greater than a threshold is not detected, the exact signal sampled 
is returned to the system. (D(Z) = 1.0) Each of the systems under 
control is handled in this fashion. When all systems have been 
sampled, the computer goes through a basic delay. This basic delay 
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is so named because it is a building block for all delays required 
for implementing variable gain amplifier solutions of any length on 
any system controlled by the computer. All required solution delays 
are formed as integer multiples of the basic delay. It is noted that 
the basic delay is much longer than the time required to sample and 
test all systems being controlled. 

When a step input is received by one of the systems, an error 
flag for the appropriate system is set and the computer exits from 
the Sense Loop to a Solution Routine. In the Solution Routine, the 
first control voltage for the flagged system is transmitted from the 
computer, and all other system flags are checked for possible simul- 
taneous step inputs. If not, the computer goes through one basic 
delay and then senses all other systems for possible step inputs. If 
no other systems have received signals, the computer returns to the 
start of the Solution Routine and begins to count the number of times 
the basic delay has been entered. After each pass through the basic 
delay, all other systems are checked for possible step inputs. When 
the pass count reaches a number which indicates that the first control 
voltage has been in the proper length of time, the computer computes 
and outputs the second control voltage, and a new pass count begins. 

If, at any time during the solution, another system 'receives a step 
input, it is flagged, and on the next pass through the Solution Routine, 
a solution similar to the one explained above is begun. 

When a variable gain amplifier solution for one of the systems 
is completed, the computer jumps to a Return Routine. In this Return 
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Routine all counters for the appropriate system are zeroed and the 
error flag is set to zero. The system is then bootstrapped by the 
procedure previously explained, and returned either to the Sense Loop 
or the Solution Routine, if the error flag of another system is set. 

In any case, after being bootstrapped, the system returns to a mode 
in which the characteristics of the control function are D(Z) = 1.0, 
and T = the basic delay. 

The characteristics of the final control philosophy may be 
summarized by stating that any single system which is being controlled 
operates in a state in which D(Z) = 1.0, and T = the basic delay when 
there is no step input above the threshold level. Thus, for small 
step inputs and noise inputs, the system would be controlled in the 
above mode. When a step input greater than the threshold level is 
received, a variable gain amplifier solution is executed and the 
control mode then returns to that specified above. In the time 
sharing mode, each system controlled is operated in this same manner, 
with the option of making the solution sampling period any integer 
multiple of the basic delay. 

To test the final control philosophy, Program 6 was written. 

For purposes of this test, two systems were time-shared. One of the 
systems used was the real system used in previous tests. The second 
system was an analog simulation of a plant having the following 

characteristics: G(S) = -77 ^ tv • The basic delay was established 

as 0.025 seconds. The solution sampling rate for the real system 
was 0.3 seconds; for the analog system, 1.0 seconds. The only 
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restriction for Program 6 was that it was written to control only 
two systems simultaneously. The basic delay time may be chosen 
arbitrarily, and the solution sampling rates may be any integer 
multiple of the basic delay. Computation of h^ and h^ was accomplished 
in the manner of Program 5. The hardware arrangement for the real 
system was the same as that employed in conjunction with Program 5, 
and the summing amplifier arrangement in and out of the converter for 
the analog system is identical to that of the real system. 

A flow chart for Program 6 is shown on the following pages. It 
is broken down into three sections: Sense Loop, Solution Routine, 

and Return Routine. A complete text of Program 6 is presented in 
Appendix VI. 

To simplify the flow chart for the Solution Routine, the negative 
step solution for System 1, and the positive step solution for System 
2 have been omitted. The omitted processes are identical in form to 
those shown. Furthermore, the "D/A h^" blocks, shown in the flow 
chart as predefined processes are handled exactly as in Program 5. 

Before presenting the flow chart for the Solution Routine, the 
following indices are defined: 

I = the number of passes through the basic delay for the 
current h for System 1. 

J = the number of h's transmitted to System 1. 

P = the required number of passes through the basic delay to 
achieve the specified variable gain amplifier solution 
sampling rate for System 1. 

K = the number of passes through the basic delay for the 
current h for System 2. 
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L = the number of h's transmitted to System 2. 

Q = the required number of passes through the basic delay to 
achieve the specified variable gain amplifier solution 
sampling rate for System 2. 



The flow chart for the Solution Routine is shown on the following 



page. 
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Program 6 
SOLUTION ROUTINE 
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The Return Routine is shown on the following page. The Return 
Routine will be shown for System 1 only. That for System 2 is identical 
except for the fact that different indices are zeroed. 
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2.3 Results 



The results of Program 6 were excellent. Both systems were controlled 
with predicted speed and accuracy. The real system response was equivalent 
to that obtained from utilization of Program 5. It is interesting to 
note that the system simulated on the analog computer exhibited exact 
theoretical response to the variable gain amplifier solution. This, of 
course, is due to the fact that there were no non-linearities in the 
simulated system. However, the exact theoretical solution obtained does 
validate the final control philosophy employed. 

The time-share program was tested for all possible time combinations 
of step inputs to the two systems. In no case was system response de- 
graded once the solution was started. In the worst case, there was a 
25 millisecond delay between a step input and the start of a solution. 

This occurred when a step was entered just as the computer started a 
pass through the basic delay in the Sense Loop. The delay, in this worst 
case, which is a random occurrence, was not discernable on the Mark II 
Brush Recorder used to measure response. 

On the following pages, the response curves for the two systems 
controlled in the time-share mode by Program 6 are shown. Also shown 
are the theoretical response curves for the real system. 

2.4 Extension of Results 

Due to the limitations imposed by the D/A, only two systems could 
be simultaneously controlled in the Digital Control Laboratory. To 
fully substantiate the validity of the final control philosophy, one 
must consider the feasibility of this philosophy in controlling more 
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than two systems. Let us consider this feasibility with respect to the 
equipment used throughout the laboratory tests, assuming only that the 
number of channels in the converter may be arbitrarily expanded. (The 
authors feel that the problems associated with more difficult tasks in 
the computer field are too often dismissed with a bland inference that a 
bigger, faster computer can do any job. The dollar economics of the 
control problem preclude this approach.) 

Using Program 6 as a basis for time-sharing control of many systems, 
the CDC-160 has sufficient memory capability to handle 12 systems 
simultaneously. By sub-routinizing wherever possible, this number could 
probably be increased to 15. For 15 systems, however, the sampling 
time, governed primarily by the analog to digital conversion time of 
120 microseconds, becomes an appreciable proportion of the basic delay 
time. In the Sense Loop, for example, the time to sample A/D, check for 
step inputs, and return D/A would be approximately two milliseconds 
compared to the 25 milliseconds of the basic delay. This added delay 
in the Solution Routine would produce significant errors in the variable 
gain amplifier solution. If the basic delay could be satisfactorily 
increased to 50 or 100 milliseconds, the two milliseconds required to 
sample the 15 systems would become insignificant and not seriously 
degrade the solution. The authors feel that employment of a 50 or 100 
millisecond basic delay, which limits the fastest solution to a step 
input to 100 or 200 milliseconds, is not unreasonable in most practical 
applications . 

The foregoing discussion has made one tacit assumption. That 
assumption is that the number of simultaneous inputs is limited to three 



52 



of the 15 systems. This is due to the 900 microsecond multiply time 
required by the CDC-168. This, of course, excludes the desired 
generality for handling any step input at any time for any system. 

With the equipment used, the 160 in conjunction with the 168, there 
is no way to compensate for this problem. 

At this point it is worthwhile to look at the advantages of using a 
different computer. Specifically, the CDC-160A would adequately solve 
the problem mentioned above. The 160A is a computer quite similar to 
the 160 in all respects with the additional capability for fast 
multiplication and division. Economically, the 160A is comparable in 
cost to a 160-168 tandem, and makes a much more compact unit. 

On the basis of the foregoing tests and theoretical extension, 
the authors feel that with a small, relatively inexpensive computer, 
such as the CDC-160A, it is quite feasible to control up to 15 systems 
simultaneously. The control capability will vary with the capabilities 
of the brand and type computer selected. The main thesis, however, 
is that an excellent, multi-system, digital control scheme may be 
implemented without going to the larger general purpose computer. 
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CHAPTER III 



3.1 Conclusion 

It was noted in the Introduction that the basic goals of this 
thesis were the controlling of hardware by digital methods, and the 
developing of techniques for implementing time-sharing control of many 
systems from the same computer. The results obtained have been detailed. 
In this section the most important findings will be summarized and their 
relevance to the general feasibility of digital control systems will 
be noted. 

The most significant aspect of controlling the hardware was the 
necessity for going to the "continuous" mode (D(Z) = 1. and a small T) 
to bring the error to zero at the conclusion of the solution to a step 
input. This characteristic, caused by the non-linearities of the system, 
would seem to be a significant disadvantage to this type control. 

However, this was not the case. In the worst case, the servo position 
was within ten percent of the desired final position at the moment of 
the switch to the "continuous" mode. Furthermore, at the switch time, 
the servo velocity was in the proper direction at a decreasing magnitude. 
This condition permitted the system to settle rapidly to the desired 
position. 

Of course, to realistically evaluate this control method, one 
must establish a performance criterion as the basis for comparison with 
other control methods. The criterion selected was one which measured 
minimum time to reach a position within ten percent of final value and 
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to remain within this ten percent boundary while settling. On the 
basis of this criterion, the time for a unity feedback continuous 
system to meet the conditions was 1.2 seconds, and 800 milliseconds 
with the addition of optimum tachometer feedback. These figures were 
obtained with the amplifier gain potentiometers set at the same levels 
used for the digital operation. By increasing the amplifier gains, a 
time of 500 milliseconds was obtained with tachometer feedback. Using 
a variable gain amplifier solution sampling rate of 300 milliseconds, 
a time of 600 milliseconds was required to meet the criterion. The 
relatively low saturation level of the amplifiers precluded the use of 
sampling rates of 100 or 200 milliseconds, but on the basis of responses 
obtained at rates of 300, 400, and 500 milliseconds, it follows that 
improved amplifiers would allow a solution time of 200 milliseconds. 

In general, the fastest speed of response obtainable by the variable 
gain amplifier method is limited by the saturation level of the 
amplifiers. Similarly, the response speed obtainable by a bang-bang 
technique is limited by this saturation level. 

On the basis of these results, the performance characteristics 
of the digital control method indicate that this method is competitive 
with other control methods. With improved amplifiers, response speeds 
from the digital methods can be much faster than those obtainable from 
continuous methods. Optimum bang-bang techniques approach the response 
speeds obtainable with the variable gain amplifier method. 

There are, of course, many criteria by which control systems 
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may be judged. The criterion above is one of speed of response. 

Although other criteria were not evaluated, it is noted that the 
digital technique is, in general, quite adaptable. For example, the 
basic hardware and software employed by the authors is readily adapt- 
able to a criterion which calls for minimum fuel expenditure. 

The results of the time-sharing program have been fully discussed 
and the feasibility of extending the program to control many systems 
has been outlined. The basic approach of utilizing a small computer 
has precluded a discussion of a further extension of the technique. 

That is, using a large, general purpose computer as the digital 
controller, with other computing tasks being time-shared with the 
control function. This is feasible with no change in the basic program 
heretofore used in conjunction with the CDC-160. In examining this 
program, one can readily see that the vast majority of time is spent 
in the basic delay block. Time spent in the basic delay is ideal for 
carrying out other computing tasks. This could readily be implemented 
on a general purpose computer with a real time clock and interrupt 
capability, and an executive routine to direct the proper sequence of 
computing tasks not associated with the control task. The ramifications 
of this large computer capability are indicated by Slaughter and 
Lackowski in a paper presented to the 1963 National Convention on 
Military Electronics in which they state: "In those cases in which a 

digital computer is available, digital control requires less hardware 
than does conventional control." M The presumption accompanying this 
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statement w as that the digital computer was required for other tasks, 
and the controlling tasks could be satisfactorily time shared. 

In view of the performance characteristics of the variable gain 
amplifier method and the time-sharing capabilities developed, the 
authors have come to the conclusion that for many applications, digital 
control methods are presently feasible and superior to other control 
methods. Due to the present cost of computers, these applications are 
presently limited to large processes or areas where the advantages of 
digital control are overwhelming. A specific example of such an 
application is in large processes such as those in the chemical or 
petroleum industries in which multiple systems must be simultaneously 
controlled. Similarly, a single large process with multiple inputs 
may be controlled by digital methods. A military application would 
be the extension of the NTDS to include the servo control of gun 
mounts and missile launchers from the master computer. A final 
application might be in the aerospace field in which the new techniques 
for optimum control on the basis of minimizing a given cost function 
are quite adaptable to digital techniques. 

Although immediate feasibility for digital control is limited 
to processes of the type mentioned above, the decreasing size and 
cost of digital computers increase their field of application. 

The investigation for this thesis was limited in scope to permit 
the extensive study of a single method. The results from this study 
may, in general, be extended to other digital control techniques. The 
authors feel that the results of this study of digital control validate 
its basic feasibility, and indicate the value of further study in the 
field. 
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APPENDIX I 



DESIGN OF DIGITAL CONTROLLERS AND SIMULATION RESULTS 



In addition to the design of a digital controller for a step input 
to a g (s+1) pl ant > a digital controller for a unit ramp input was also 
designed. Referring to Figure 1 and Figure 2 in Section 1.2 we arrive 
at the following equations. 

T = 1.0 

X 1 £(k+l)Tj = |^~0-368K^J x^(kT) + 0.632x2 (kT) + (0. 368K^)r(kT) 
x 2 [ ( k + l)T] = (-0.632K k )x 1 (kT) + 0.368x 2 (kT) + (0.632K k )r(kT) 



For a unit ramp input r(kT) = k. 



k = 0; x^O) = x 2 (0) = 0.0 
(T) = 0 ; x 2 (T) = 0 

k = 1 

x^ (2T) = 0.368^ 
x 2 (2T) = 0.632^ 

To obtain a solution we must set x^(nT)=r(nT) and x 2 (nT)=l. We there- 
fore must proceed to the next period in order to solve for the given 
conditions . 

Therefore k = 2 

x 1 (3T)= |^l-0. 368K 2 J (0.368^) + (0.632) 2 K 1 + 2(0.368K 2 > = 3.0 
x 2 (3T)= f-0.632K 2 l (0.368^) + (0.368) (0.632^)4-2 (0.632K 2 ) = 1.0 
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from which 



K x = 3.82 ; K 2 = 0.31 

and 

h(T) = 3.82 ; h(2T) = 0.183 



It should be noted that for this system to maintain a zero error after 
it has arrived "home" a constant output of h(nT)=l must be fed into the 
continuous system. This will keep the system moving at a unit velocity 
and thus maintain a zero system error. 



for a 



A digital controller using the same design theory was developed 
1 



plant in response to both unit step and ramp inputs. 



(s+l) (s+2) 

The derivations and explanations follow. The procedure used is the 
same as that in Section 1.2. 

Plant: G(S) = (s+1 ) ( s+2 ) 

Sampling Period: T = 1.0 

Block Diagram: 








1 






D(S) 


1 

A ! 

T | 
1 


z . o.h. 


h III 











(s+l) (s+2) 



Xj^(t) 



D(Z) 



Figure 24 
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-2 




^ j^(k+l)Tj = J^2 e T -e 2T -I^(.5-e T +.5e 2T )J x ]L (kT)+(e T -e" 2T )x 2 (kT)+ 



^(.S-e - ^ . 5e” 2T ) r(kT) 



[(k+l)Tj = jle'^-Ze'^Ce^-e' 21 ) Jx x 



e A -e i " L ) | x n (kT)+(2e T -e T )x 2 (kT)+ 



-T -2T 

\( e -e Zi ) r(kT) 



For T = 1.0 



x 1 |^(k+l)Tj = |T 601- . 19951^^ x^kT) + .233x 2 (kT) + .19951^ r(kT) 
x 2 £(k+l)T J = jj- .466 - . 233kjj x^kT) - .098x 2 (kT) + .2331^ r(kT) 



(I) Let r (kT) =1.0 

k=0 ; x^CO) = x 2 (0) = 0 

x n (T) = . 1995K ; x„ (T) = .233K 

1 O 2 o 

These two equations clearly cannot be solved for the final conditions 
of: x 1 (T) = 1.0 ; x 2 (T) = 0.0 

Therefore let k = 1 

x, (2T) = 0.1742K - 0.0398K K n + 0.1995K, = 1.0 

1 o o 1 1 

x 0 (2T) = -0.11575K - 0.0465K K n + 0.233K, = 0.0 

2 o o 1 1 

from which: 

K = 3*65 ; K, = 6.73 
o 1 

and 

h (0) = 3.65 ; h(T) = 1.832 
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Note: Since the plant involved is of type zero it will have an inherent 



steady-state error to a step input. To eliminate this error a "steady 
state h" will have to be put into the system after the error has been 
initially reduced to zero. To determine this "steady state h" the final 
value theorem was used as follows: 

Assume: h(t) = at 4- b 

H(S) --!+-! 



x 1 (s) = 



bs+a 



<s+l)(s+2) 



x l (c ° > dl 



(bs+a)e 

2 

s -f 3s4-2 



st 



s = 0 



x iM = 



2b+2at-3a a ^ ^ 



or 



= 0 ; b = 2 



/. h (nT) =2 n « 2,3. 



(II) Now let r(kT)=k i.'.e. unit ramp 

In this case the final values desired are: x^(nT)=k; X2(nT)=1.0 
for n the smallest possible integer. 

For: k = 0; x^ (0)=X2 (0)=0 

x^D-o ; x 2 (t) = o 

Therefore let k=l 

x 1 (2T)= 0.1995KL 1 

x 2 ( 2 T )= 0 . 233<. 1 
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These two equations cannot be solved so we must proceed to fe=2. 



x 1 (3T)= 0.17425^ - 0.0398^^ + 0.399K 2 = 3.0 

x 2 (3T)= -0.11575^ - 0.0465K 1 K 2 + 0.466K 2 = 1.0 
From which 



= 7.85 ; K 2 = 18.71 

and 

h(T) = 7.85 ; h(2T) = 8.14 



Again to keep the output equal to the input for this type zero system 
we must output a "steady state h". From the final value theorem this 
is : 



h(nT) = 2n+3 n 3,4,5,.... 



The following pages contain the graphical simulation results for 

1 



the type one system (G(S) = 



s (s+1) 

responses of a type zero system (G(S) = 
and unit ramp inputs. 



) to a unit ramp input and the 

1 



(s+l) (s+2) 



) to both unit step 



63 



OGQ 001 OQZ 03Q 00 l i OQZ OQG 007 




— WM. E - J.00E4W LMJTS/JNCH. 

'1-.5MUE - ]ML-*tiZ LNITVJNCK 

BROWNE DIGIT ,1R G < s > - ustt r < kT > - k 
RUN 1 OUTPUT US. TIME 

Figure 26 
64 




?cer oca 32^0 cz0 3^0 ota 300 000 



\ 








































\ 












\ 












\ 






• 












i 

i 


L 













ma J 0 J 3 Q aaj 004 aos 



*-SMLIE - j^U 3 £-+J 0 ja UMJTVJMCIi. 
M-.SMLK - 5 J& 0 IE-J 8 J UNITS/JNCH 



BROUIME DIGIT 1R 

RUfH I UELOCITY US.TIME 



Figure 27 
65 



-0 10 -00£ 000 0OZ On 0 




*-3MLIt - UTUT5/JNCM 

'l-SCfllR - 5J&0E-A.1 UNJT-5/JNCH 

BROWNE DIGIT 1R 

RUN J EDOT US. ERROR 

Figure 28 
66 



000 OO'Z OO'i OOG 00c) 010 Oil 011 




K-3MLE - JJWJE-+J&B UNJTWJNCM 

'fOMLE - Q.&QE-J&J UNJTvVJNCH. l 

BROWNE DIGIT 2 G(s) = ( ^nt^ n r ( kT >= i -° 
RUN 1 OUTPUT US, TIME 

Figure 29 
67 



Tii o 3i o ozo eco eoo wo wo ooo 




K-:CFLZ - UNJTS'JNCh. 

'f-JCflLE - Q.&0E-0J LJNJT5/JHCH. 



BROWNE DIGIT 2 

RUN 1 UELOCITY US.TIME 



Figure 30 
68 



-010 -003 -005 -001 -002 000 002 001 




X--5MLE - O.M£-ei UNJTS/JNCH. 
'l-dCPLE - UKIT5/JNCH. 



BROUIME DIGIT 2 

RUN 1 EDOT US. ERROR 



Figure 31 
69 



000 002 000. 003 004 0OZ 00G 00? 




hm mi mi m3 m\ ms 



x-$mlic « j^-+aei units/ jmch. 

't-SMUE ' LMIT3/ JriCM 

BROklME DIGIT 2 R a(s) - wlw *■<«> - k 
Rm 1 OUTPUT US„ TINE 

Figure 32 
70 



000 80S, 020 015 02. 0 025 030 835 




x-scalz - j.&2£4jaa ukus/jmcm. 

^-5Ca£ - 5jae£-£ij UTHirS/JNCn. 

BROWNE DIGIT 2R 

RUM 1 UELOCITY US.TIME 

Figure 33 
71 



~01Q ~0@£ 000 00S 010 0 1& 




H-scpLfc - 'i.w-aj units/incm 

'l-.SMLiE « 5J&0£-£)J UNJTS/JNCtt 

BROWNE DIGIT OR 

RUN 1 EDOT US. ERROR 

Figure 34 
72 



APPENDIX II 



SIMULATION PROGRAM 

The following pages contain the simulation program used in 
verifying the digital controller design for a "deadbeat response" to a 
unit step input. The programs for simulating a unit ramp response are 
similar. 
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..J03127F, BROWNE, E.R. 

PROGRAM DIGIT 1 

DIMENSION X ( 30 ) ,XDOT (30) »C( 15) 
C ( 1 0 ) = 1 . 0 

1 CALL INTEG1 (T,X,XDOT,C) 

INPUT =1.0 

IF { T- 1 « *C ( 1 ) ) 10,11,12 

10 HOLD® C(2) 

GO TO 14 

11 HOLD = C ( 3 ) 

GO TO 14 

12 IF (T-2.*C(1>) 11,13,13 

13 HOLD = C ( 4 ) 

14 XDOT ( 2 ) =HOLD-X ( 2 ) 

XDOT ( 1 ) =X < 2 ) 

ERROR= 1 • 0-X ( 1 ) 

EDOT=-X ( 2 ) 

X ( 3 ) =ERROR 
X ( 4 ) = E D 0 T 
C( 11 ) =50 • 0 
GO TO 1 
END 
END 
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BROWNE DIGIT 1 

ONE RUN IS CALLED FOR 



INPUT DATA RECORD 



ORDER OF EQUATIONS = 2 

INITIAL TIME * . OOOO i + CO 

FINAL TIME = . SOOOE+CC 

STEP SIZE = .2000E-C3 



THE NCN-ZERO CONSTANTS 



C( 

C( 

C( 



1 ) = 

2 ) = 

3) = - 



C(I), ARE 
1 OOOE+CO 
1 056E+C3 
9580E+C2 



ALL THE INITIAL CONDITIONS ARE ZERO 



THE COLUMN HEADINGS AND THE CORRESPONDING VARIABLES ARE, 



TIME 


X{ 


0) 


OUTPUT 


X{ 


1 ) 


VELOCITY 


X( 


2) 


ERROR 


X( 


3) 



THE INDIVIDUAL GRAPH TITLES AND THE CORRESPONDING VARIABLES ARE 



OUTPUT VS. TIME 


X( 


I ) 


VS. X< 


0) 


VELOCITY VS. TIME 


X( 


2) 


VS. X( 


C) 


EDOT VS. ERROR 


X( 


4 ) 


VS. X( 


3) 
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BROWNE DIGIT 1 



TIME 



• OOOOOE+OO 
. 10000E-0 1 

• 20000 E-0 1 
. 3Q0QOE-O 1 
. 40000 E-0 1 
.50000 E-0 1 
. 60C00E-0 1 
. 70000 E-01 
. 80000 E-0 1 
. 90000 E-0 I 

. 1 0000 E+00 
. 1 1000E+00 
. 12000E+00 
. 13000E+00 
. 14000E+00 
. 15000E+00 
. 16000E+00 
. 17000E+00 
. 1800CE+00 
. 19000E+00 

. 20000 E + 00 
. 2 1 000E+00 
. 22000E+00 
. 23000E+00 
. 24000E+00 
. 25000E+00 
. 26000 E+00 

* 27000E+00 
. 28000E+00 
. 29000E+00 

. 30000 E + 00 
. 31 OOOE+OG 
. 32000 E+00 
. 33000 E+00 
. 34000 E + 00 
. 35000 E + 00 
. 36G00E+00 
. 37000E+00 
. 38000 E + 0 0 
. 39000 E + 00 

.40000 E+00 
. 4 1000E+00 
.42000E+00 
.43000 E+00 
.44000 E+00 
. 45000 E+0 0 
. 46000E+00 
.47000E+00 
.48000 E+00 
. 49000E+00 



OUTPUT 



. 00000 E + 00 
.52624E-02 
. 2C980E-0 1 
.470 4 8 P-01 
.83365E-C 1 
. 12983E + 00 
. 1 86 33E+00 
. 25279E+00 
. 32909E+00 
.4151 3E+00 

.51083E+00 
. 60598E+00 
.69065E+00 
. 76495E +00 
•82897E+0C 
. 88283E+00 
.92662E+00 
. 9 60 4 4E +00 
.98439E+00 
. 99857E+C0 

. 10031E+01 
. 10027E+C1 
. 10024E+G1 
. 10021E+C1 
. 10018E+01 
• 100 1 56 + 01 
. 10012E+01 
. 10008E+01 
. 1C005E+01 
. 1 0002E + 0 1 

. 99993E+00 
.99964E+00 
. 99934 E+CO 
. 99905E+00 
. 99876E+00 
.99848E+00 
.99820E+C0 
. 9979 2E + 00 
. 99764E+00 
. 99737E+C0 

. 997C9E + C0 
.9968 3E +00 
.99656E+C0 
. 99630E+00 
. 99604E+00 
.99578E +00 
. 99552E + 00 
. 99527E+0C 
. 99502E+00 
.99477E+00 



VELOCITY 



.00000 E+CO 
. 10507 E + 01 
.209 10 E + 0 1 
.312106+01 
. 4 1 406 E + 0 1 
.5 1 502 E + 01 
.614 97 E+01 
.71392E+01 
.81 1 89 E + 0 1 
. 90889 E+01 

. 1 C042E+C2 
.898 93 E+01 
.79466 E + 01 
.691 43 E+01 
. 58923E + 0 1 
.48805 E+01 
.38787 E + 01 
.288 69 E + 01 
. 1 9049E+C1 
.93272 E+CO 

- . 29786E-01 
-.32651 E-01 
-.32326 E-01 
-.32005 E-01 
-.3 3686 E-01 
-.31371 E-01 
-.3 1 059 E-01 
-.30750 E-01 

- . 3G4 44 E-0 1 

-.30141 E-01 

-.29841 E-0 1 
-. 29544E-0 1 
-. 29250 E-01 

- .28959 E-01 
-.28671 E-01 
-.28386E-01 
-. 28 1 03 E- 0 1 
-.27824 E-01 
-.275 47 E-01 
-. 27273E-01 

-.27001 E-01 
-.267 33 E-01 
-.2 64 67 E-01 
-.26203 E-01 
-. 25942 E-0 1 
— . 2 5684 E-0 1 
-.25429E-01 
-.25176 E-01 
-.24925 E-01 

- .24677E-01 



ERROR 



. 1 OOOOE + G 1 
. 99474E+00 
. 979C 2E + 00 
. 952956+00 
. 91 664E+00 
. 870 1 7 E+00 
. 81367E+00 
. 74721 E+00 
. 6709 1E+00 
. 584 87 E+00 

. 4891 7E+00 
. 39402E+00 
.30935 E+00 
. 235C5E+0C 
. 17103E+C0 
. 1 1 717E+00 
. 73384E-01 
. 39565E-01 
. 1561 4E-01 
. 14340E-02 

-. 30727E-02 
- . 2744 5E-02 
-.241 96E-02 
- . 2098CE-02 
-. 1 77 95E-02 
-. 14643E-02 
-. 1 1521E-02 
- . 843C7E-03 
-.5371 0E-03 
- . 234 1 8E-03 

. 65725E-04 
. 362 65E-03 
. 65662E-03 
. 94766 E-03 
. 1 2358E-02 
. 1521 IE-02 
. 1 8035E-02 
.20832E-02 
. 236CCE-02 
. 2 634 1 E-02 

. 29055E-02 
.31741 E-02 
. 3440 1 E-02 
• 37035E- 02 
. 39642E-02 
. 42223E-02 
. 44779E-02 
. 473 09E-02 
.498 14 E-02 
. 52294E-02 
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BROWNE DIGIT 1 



TIME OUTPUT VELOCITY ERROR 

. 50000 E+OO . 99453E + 00 -.24432E-01 . 54750E-02 



NORMAL 


STOP AT 


FINAL TIME 










GRAPH 


TITLED . 


• 


BROWNE 


DIGIT 


1 












RUN 1 






OUTPUT VS 


. TIME 


GRAPH 


TITLED . 


• 


BROWNE 


DIGIT 


] 












RUN 1 






VELOCITY 


VS. TIME 


GRAPH 


TITLED . 


• 


BROWNE 


DIGIT 


1 












RUN 1 






EOOT VS. 


ERROR 



THE ONE RUN CALLED FOR HAS BEEN COMPLETED. 



STOP 

TIME, 1 MINUTES AND 14 SECONDS 
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APPENDIX III 



CONTINUOUS PLANT 

The continuous plant which was to be controlled in a sampled- 
data manner consisted of two operational a.c. amplifiers, a d.c. power 
amplifier, and a d.c. torque motor manufactured by the Inland Motor 
Corporation of Pearl River, New York. 



h (t\ 



a.c. 



Amplifiers 



Figure 35 

As is seen in the figure, the control signal to the plant was 
received from the zero-order hold and fed into the first operational 
amplifier. The second operational amplifier transmitted the signal to 
the d.c. power amplifier where the signal size was increased further. 
From the d.c. power amplifier the control signal was fed to the armature 
of the motor. 

The Inland motor has a permanent magnet d.c. field and is armature 
controlled. An inertia disc was attached to the shaft to slow down the 
response of the system. Also on the motor shaft were a tachometer for 
velocity pick-off and a potentiometer which was used for position feed- 
back. The motor was rated at 0.47 ampere and 90 volts. 5 
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To provide the needed power to the motor an amplifier whose 
capabilities were greater than the common commercial type was needed. 

The d.c. amplifier used was the same one that had been constructed for 
use on a previous thesis using the same motor. ^5 J As seen in Figure 
35 the input to the d.c. amplifier was double ended. The amplifier 
operation is push-pull which reduces third harmonic distortion and 
provides a balanced output. (See Figure 36 for a circuit diagram of 
the d.c. power amplifier.) 

To insure that the d.c. power amplifier was operating correctly 
two tests were performed. The first test ‘checked the frequency response 
characteristics of the a.c. and d.c. combination. The second test 
plotted the output vs. the input of the d.c. amplifier alone thus 
checking the saturation characteristic of the amplifier. Test data 
and curves are shown on pages 82-84 . It should be mentioned that some 
difficulty was encountered with the d.c. amplifier. The amplifier needs 
a supply voltage of 300 volts with a corresponding current of one 
ampere. Since this power was not available directly in the room where 
the equipment was set up, it was decided to place four portable power 
supplies in parallel. This arrangement proved fairly satisfactory with 
the exception that the power supplied was not regulated. This caused 
difficulty in balancing the amplifiers and once they were balanced they 
did not remain so for any length of time. This balancing problem was 
bothersome but in no way invalidated the results or conclusions of 
this report. 

To determine the transfer function of the amplifier motor 
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combination a standard frequency response test was conducted on the 
closed loop system. The input sinusoid was supplied by a Hewlet 
Packard Low-Frequency Generator Model 202A. Unity feedback was em- 
ployed and a Mark II Brush Recorder was used to measure the input and 
output. (See Figure 37 ) 




Figure 37 

The recorded closed loop data was first plotted as a Bode 
diagram. (Magnitude and Phase) The closed loop data was then con- 
verted to an open loop transfer function by using a Nichols Chart and 

a second Bode plot. The open loop transfer function was determined to 
A 

be: G(S) = — — rr with A set at 40. (See pages 85-89 for test 

S V. S T J J 

data, Bode diagrams, and Nichols Chart.) 
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Frequency Response of Amplifiers in Cascade 



frequency 
(cps ) 


u> 

(radians ) 


input 
(volts ) 


output 
(volts ) 


gain 


0.1 


0.628 


5.4 


55.0 


10.2 


0.2 


1.255 


5.4 


55.0 


10.2 


0.3 


1.885 


5.3 


55.0 


1C. 4 


0.4 


2.510 


5.3 


55.0 


1C. 4 


0.5 


3. 140 


5.2 


55.0 


10.6 


0.6 


3.770 


5.2 


55.0 


10.6 


0.7 


4.390 


5.2 


55.0 


1C. 6 


0.8 


5.020 


5.2 


55.0 


10. 6 


0.9 


5 . 650 


5.2 


55.0 


10.6 


1.0 


6.280 


5.2 


55.0 


10.6 


1.1 


6.900 


5.2 


55.0 


10.6 


1.2 


7.550 


5.2 


55.0 


10.6 


2.0 


12.550 


5.2 


55.0 


10.6 


4.0 


25.100 


5.2 


55.0 


10.6 


6.0 


37.700 


5.1 


55.0 


10.8 


e.o 


50.200 


5.1 


55,0 


10.8 


10.0 


62.800 


5.0 


54.0 


10.8 


12.0 


75.500 


5.0 


53.0 


10.6 


15.0 


94.500 


5.0 


53.0 


10.6 


20.0 


125.500 


5.0 


52.0 


10.4 


40.0 


251. COO 


4.6 


48.0 


10.4 


50. 0 


314.000 


4.6 


42.0 


9.2 


70.0 


439.000 


4.2 


34.0 


. 6.6 



Saturation Characteristic of d.c. Amplifier 
frequency's 0.7 cps; u) - 4.39 



input 
(volts ) 


output 
(volts ) 


0.3 


1.0 


1.3 


4.4 


3.1 


10.0 


5.6 


17 c 5 


7.6 


24.0 


10.0 


33.0 


13.0 


43 . 0 


15.5 


50.0 


18.0 


58.0 


20.0 


64. C 


22.0 


71.0 


25.0 


77.0 


27.0 


82.5 


30.0 


85.0 


32.0 


90.0 


34.5 


95. ‘ 0 


38.0 


100.0 


42.0 


102.5 


40.0 


J 05 0 


so.o 


1 07 - 5 " 


52 . 0 


l 1 0.0 




i.o to.o ioo.o 

oo CRADiANS/jecoND') 



Frequency Response Data 



frequency 

(cps) 

\ 


CcD 

(radians ) 


magnitude 
(db. ) 


phase angle- 
(degrees ) 


0.1 


0.628 


0.0 


- 0 I 


0.2 


1.25 


0.7 


- 11 


0.25 


1. 57 


1.3 


- 16 


C. 3 


1.88 


1.6 


- 11 


0.35 


i. 2.20 


2.6 


- 22 


0.4 


2.55 


2.4 


- 25 


0.45 


2 . 62 


3.18 


-■ 23 


0.5 


3.14 


3.5 


- 33 


0.55 


3.45 


3.7 


- 29 


0.6 


3.78 


3.75 


- 31 


0. 625 


3.93 


3.75 


- 37 


0.65 


4. 08 


4.3 


- 37 


0.675 


4.25 


3.75 


- 39 


0.7 


4.40 


5.2 


- 41 


0.725 


4.56 


5.2 


- 37 


0.75 


4.71 


5.2 


- 43 


C. 775 


4.87 


3.75 


- 45 


0.6 


5.03 


4.3 


- 47 


0.825 


5.19 


5.75 


- 48 


0.85 


5.34 


6.85 


- 86 


0.675 


5. 50 


6.7 


- 89 


0.9 


5.65 


6.3 


- 92 


0.95 


5.96 


5.3 


-102 


1.0 


6.28 


4.25 


-106 


1.05 


6.59 


3.25 


-112 


1.1 


6.91 


2,6 


-117 


1.5 


9.45 


-2.9 


-124 
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APPENDIX IV 



MEASUREMENT OF QUANTIZATION ERRORS 

This appendix will cover the determination of quantization errors 
in the sampled-data control system. It will only deal with the real 
plant since the main thesis effort was directed to the control of this 
plant. The theory used was described by John B. Slaughter of the Navy 
Electronics Laboratory in the January 1964 IEEE Proceedings on Automatic 
Control. |^6 

Quantization, or round off, errors occur because of the inevitable 
conversion operations required when analog and digital devices are 
connected together in a sampled-data closed loop system. Quantization 
is the process of converting a signal in analog form to its digital 
approximation. This is best illustrated by Figure 44. 




The dashed line in the figure is the desired linear response, while 
the staircase represents the actual output of the quantizer. The 
range of input magnitudes is divided into disjoint intervals (h^) which 
are not necessarily equal. All magnitudes falling within an interval 
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are equated to a single value within the interval. This results in 
the digital approximation to the analog input. It can be seen that the 
maximum round off error of the quantization operation is h/2. (h is 
called the quantizing level.) 

Quantization is non-linear in the sense that a discrete set of 
amplitude levels is produced for some continuous range of input. Be- 
cause of this an analytical expression for the effects of this error 
is difficult to obtain. However, by assuming that the quantizing 
levels are kept small the concepts of probability may be used. The 
quantizer is treated as a summing point with uniformly distributed 
noise introduced at that point. 

In continuing, we will refer to the real plant which was used in 
the sampled-data control problem. 




In Figure 45 we have the system that was used in this study. The 



sampling period is 0.3 seconds, and the digital controller is designed 
to give a ’’deadbeat response" to a step input. From z-transform theory 



we have: 



G(Z) 

D(Z) 



.368Z~ 1 +.264Z~ 2 
1 - 1.368Z' 1 +.368Z 
1.582 - . 582Z' 1 
1 + .418Z" 1 



( 1 ) 

( 2 ) 
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The equations describing the system may also be written as: 



x(k + 1) = Ax(k) + Dr(k) (3) 

c (k) = Bx(k) (4) 

Now if we decompose our system into its canonical form equations 
2 and 3 may be written by inspection. 




Canonical form of sampled-data control system 



Thus, if we define our states as shown in Figure 46 we arrive at: 



A = 



D = 



0 

-.831 
■ 1.01 
"o “ 

.42 



1 

.833 

-1.365 

( 6 ) 



0 

-.349 

-.424 



(5) 



[ 



B = 1.01 1.365 0 



] (7) 



Before the input, r(k), or the output, c(k), are introduced into 
the digital controller they pass through the A/D converter and are 
changed from analog signals to digital signals. Also, since the 



92 





digital computer must operate with finite word lengths 5 round off occurs 
for each word the computer processes or outputs. Therefore quantization 
operations occur and should be accounted for wherever these processes 
happen. This may be done by introducing quart izers in the canonical 
form of the system as shown in Figure 47 . 




Figure 47 



We now introduce a new vector R(k) which is the difference between 
the quantized and un-quantized terms. The quantized system may now be 
described by: 

x (k+l)= Ax (k) + R(k) + Dr(k) (8) 

-q “q - 

where the subscript q indicates the quantized system. It may be 
seen by inspecting the above equation that if the un-quantized system 
is stable then the quantized system is also stable. 

To obtain the difference between the quantized and un-quantized 
systems we define: 

e(k) = x q (k) - x(k) (9) 
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from which: 



( 10 ) 

( 11 ) 



£(k+l) = Ae_(k) + R(k) 
c q (k) - c(k) = Be(k) 

Thus in our problem: 

e 1 (k+l) = e 2 (k) (12) 

e 2 (k+l)= -.831e (k) + .833e 2 (k) - .349e 3 (k) + R (k) (13) 

e (k+l)= -l.Ole (k) - 1.365e (k) - . 424e 3 (k)+R 3 (k) (14) 

We must now find an upper bound on R and hence on e. From this we 
may obtain the maximum output error due to quantization. Since R^(k) 
is the difference between the state vector component, x^(k+l), obtained 
in the quantized and un-quantized systems, an upper bound on R^ (k) is 
the greatest error due to quantizing in one sampling period. 

In our system we have: 



R x (k) = 0 



R 2 (k) £ 1*42^ + .21h 2 + h 2 h 3 | = H 2 



R 3 (k) = 



h, + .5h„ 1= H 



(15) 

(16) 
(17) 



where h^/2 is the A/D quantizing error, h 2 /2 is the computer 
round off error, and h 3 /2 is the D/A quantizing error 
To find an upper bound on e(k) we must convert equations 12, 13, and 
14 in their z-transform equivalents. 

Thus : 

e x (z) = z" 1 e 2 (z) • (18) 
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(19) 



e^(z)= -.831z ^e^(z) + .833z ^e^(z) 

e^(z)= -l.Olz ^e^(z) - 1.365z ^e^Cz) 

Referring to equations 7 and 11 it is seen that to find the overall 

system error due to quantization we only have to solve for e^(z) and e^(z) . 

Also, the steady state error e^ equals the steady state error e^ . That 

is: lim (1-z '*‘)e 1 (z) = lim (1-z ^)e (z) 

z*l 1 z+1 1 

Therefore after solving equations 18, 19 and 20 simultaneously and 

applying the final value theorem we arrive at: 

1.424H - .349H 

e (oo) = e (oo) = * (21) 

.595 

and since in the hardware used h^ = 2.44 millivolts, 
h^ = .0244 millivolts, and h^ = 2.44 millivolts 
then 

(oo ) = e^ (oo) = 15.3 millivolts (22) 

Therefore c (k) ~ c(k) = 36.35 millivolts 

This difference between the quantized and unquantized systems was less 
than 1% of the final value of the system for a five volt step which was 
ordinarily used. Therefore no problems were encountered with quantization 
errors in the control of the real system. 



-i “a 2 ’ 1 

349z e~(z) H — 

1-z" 1 



-1 H 3 Z 

• 424z e (z)+ r~ 

1-z" 



( 20 ) 
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APPENDIX V 



PROGRAM FLOWCHART SYMBOLS 



Symbol 




□ 

o 



Represents 

Processing 



External Frsnctlon 
Terminal or Connector 



c 



) 




Decision 



Indexing Operation 



Predefined Process 



Jump 
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APPENDIX VI 



PROGRAMS FOR THE CDC-l6o COMPUTER 



PROGRAM 1 



Location 


Contents 


Comments 


1000 


7500 


Call A/D 


1001 


l 4 oi 




1002 


7600 


Read into A 


1003 


4000 


Store in 0000 


1004 


7500 


Call D/A 


1005 


2401 




1006 


7304 


Output ( 0000 ) 


1007 


0001 




1010 


6510 


Jump to 1000 


1011 


64 li 




1012 


0000 
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Program 2 



Location 


Contents 


Comments 


1000 


7500 


Can A/D 

* 


1001 


1401 




1002 


7600 


Read into A 


1003 


4062 


Store 'New value' 
in 0062 


1004 


3460 


Subtract 'Old 
Value' 


1005 


4065 


Store 'Diff' in 

OO65 


1006 


3063 


Add Increment 


1007 


6315 


If sum negative, 
jump to 1024 


1010 


2065 


Load 'Diff' into A 


ion 


3463 


Subtract Increment 


1012 


62n 


If difference positive 
jump to 1023 


1013 


2062 


Load 'New Value' 


1014 


4o6o 


Store as 'Old Value' 
in 0060 


1015 


7500 


Can D/A 


1016 


2l»01 




1017 


7306 


Output (0060) 


1020 


0061 




1021 


6521 


Jimp to 1000 


1022 


6422 




1023 


7071 


Jump to Negative Step 
Solution (0071) 


1024 


7072 


Jump to Positive Step 
Solution (0072) 


1025 


0060 
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Location 

0060 

0062 

OO 63 

OO 65 

0071 

0072 



PROGRAM 2 (CONTINUED) 
LOW CORE ALLOCATION 
Contents 



0010 

7700 

7700 



Comments 
'Old Value' 
'New Value' 
Increment 
•Diff' 

Dummy start of 
Negative Step 
Solution 

Dummy start of 
Positive Step 
Solution 
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PROGRAM 3 



Location 


Contents 


Comments 


1030 


7500 


Call D/A 


1031 


2401 




1032 


7317 


Output (0020) 


1033 


0021 




1034 


2050 


Zeroise counter 


1035 


3450 




1036 


4050 




1037 


5450 


Add one to 
counter 


io4o 


2050 


Subtract counter 
total from 


io4i 


3451 


preset total 


1042 


6002 


If difference -- 0 
jump to 1044 


1014-3 


6504 


If difference not 
= 0, jump to 1037 


1044 


7300 


Call D/A 


1045 


2401 


1046 


7304 


Output (0022) 


1047 


0023 




1050 


7700 


Halt 


1051 


0020 




1052 


0022 
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PROGRAM 3 (CONTINUED) 



LOW CORE ALLOCATION 



Location 


Contents 


Comments 


0020 


0631 


Converted as 
-4 volts 


0022 


71*16 


Converted as 
-6 volts 


0050 


— 


Counter 


0051 


3U03 


Preset total for 
0.1 second delay 



if 



101 



1001 

1002 

1005 

1004 

1005 

loo6 

1007 

1010 

ion 

1012 

1015 

1014 

1015 

1016 

1017 

1020 

1021 

1022 

1025 



PROGRAM 4 



Contents 

7500 

l4oi 

7600 

4o62 

546o 

4065 

5065 

6515 

2065 

5465 

62n 

2062 

4o6o 

7500 

2401 

7506 

0061 

6521 

6422 

7071 

7072 
0060 



Comments 



Sense Loop 
(See PROGRAM 2) 
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PROGRAM 4 (CONTINUED) 



Location 


Contents 


Comments 

NEGATIVE STEP SOLUTION 


1026 


2060 


Load 'Old Value' 


1027 


3014 


Add H° 


1050 


4020 


Store adjusted H° 
in 0020 


1031 


2060 


Load 'Old Value' 


1032 


3015 


Add H 1 


1033 


4022 


Store adjusted H 1 
in 0022 


1034 


7500 


Call D/A 


1035 


2401 




1036 


7327 


Output (0020) 


1037 


0021 




io4o 


2052 


Start Delay 


104i 


3452 


(See PROGRAM 3 ) 


1042 


4052 




1043 


2050 




1044 


3450 




1045 


4050 




1046 


5450 




1047 


2050 




1050 


3451 




1051 


6002 




1052 


6504 
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1054 

1055 

1056 

1057 

lo6o 

1061 

1062 

1063 

1064 

1065 

1066 

1067 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 



PROGRAM 4 (CONTINUED) 



Contents Comments 



5452 




2052 




3453 




6002 




6514 




7500 


Call D/A 


2401 




7304 


Output (0022) 


0023 




7073 


Jump to 1130 


0020 




0022 


POSITIVE STEP SOLUTION 


2060 


Load ‘Old. Value' 


3016 


Add H° 


4024 


Store adjusted H° 
in 0024 


2060 


Load 'Old Value' 


3017 


Add H 1 


4026 


Store adjusted H 1 
in 0026 


7500 


Call D/A 


2401 




7327 


Output (0024) 



D4 



PROGRAM 4 (CONTINUED) 



Location 


Contents 


Comments 


1100 


0025 




1101 


2052 


Start Delay 


1102 


3452 


(See PROGRAM 3) 


1103 


4052 




1104 


2050 




1105 


3450 




1106 


4050 




1107 


5450 




1110 


2050 




1111 


3451 




1112 


6002 




1113 


6504 




1114 


5452 




1115 


2052 




1116 


3453 




1117 


6002 




1120 


6514 




1121 


7500 


Call D/A 


1122 


2401 




1123 


7304 


Output (0026) 


1124 


0027 
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PROGRAM 4 (CONTINUED) 



Location 


Contents 


1125 


101 b 


1126 


002b 


1127 


0026 


1130 


2052 


1131 


3*452 


1132 


*4052 


1133 


2050 


113 *+ 


3*450 


1135 


*4050 


1136 


5*450 


1137 


2050 



11*40 


3*451 


11*41 


6002 


11*42 


650*4 


11*43 


5*452 


11*41 


2052 


11*45 


3*453 


11*46 


6002 


11*47 


651*4 


1150 


7500 



Comments 
Jump to 1130 



Start Delay 
(See PROGRAM 3 ) 



BOOTSTRAP 



106 



1153 

115*+ 

0014 

0015 

ooi6 

0017 

0020 

0022 

0024 

0026 

0051 

0053 

0060 

0062 

0063 

0065 

0070 

0071 

0072 

0073“0( 



PROGRAM 4 (CONTINUED) 



Contents 

1401 

7600 

4o6o 

7070 

LOW CORE ALLOCATION 



3403 

0003 



1000 

1026 

1067 

1130 



Comments 

Read into A 

Store as 'Old 
Value' 

Jump to 1000 

H° (-) 

H 1 (-) 

H° (+) 

H 1 (+) 

Adjusted H° (-) 

Adjusted H 1 (-) 

Adjusted H° (+) 

Adjusted H 1 (+) 

Preset total for 
0.1 second delay 
Number of passes 
through 0.1 second 
delay 

'Old Value' 

'New Value' 

Increment 

'Diff 
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PROGRAM 5 



location 


Contents 


Comments 


1000 


7500 


Sense Loop 


1001 


1401 




1002 


7600 




1003 


4062 




1004 


3460 




1005 


4065 




1006 


3063 




1007 


6315 




1010 


2065 




1011 


3463 




1012 


6211 




1013 


2062 




1014 


4o6o 




1015 


7500 




1016 


2401 




1017 


7306 




1020 


0061 




1021 


6521 




1022 


6422 


Jump to Negative 


1023 


7011 


Step Solution 
Jump to Positive 


1024 


7012 


Step Solution 


1025 


0060 
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PROGRAM 5 (CONTINUED) 



Location 


Contents 


Comments 


1026 


7500 


Negative Step 
Solution 


1027 


l4oi 


Call A/D 


1050 


7600 


Read into A 


1031 


3460 


Subtract *01d 
Value' 


1032 


4024 


Store E° 


1033 


0101 


Multiply E° 
by 100 K° 


1034 


7000 




1035 


0031 




1036 


0020 




1037 


0024 




104o 


0030 




104l 


0101 


Divide product 
by 100 


1042 


7000 




1043 


004l 




1044 


0030 




1045 


0002 




1046 


0030 




1047 


2030 


Load H° 


1050 


3060 


Add 'Old Value' 


1051 


h 030 


Store aid jus ted H° 
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PROGRAM 5 (CONTINUED) 



Location 


Contents 


Comments 


1052 


7500 


Call D/A 


1053 


2401 




1054 


7353 


Output adjusted 
H° 


1055 


0031 




1056 


2052 


Delay 


1057 


3452 





1060 


4052 


106l 


2050 


1062 


3450 


IO63 


4050 


1064 


5450 


1065 


2050 


1066 


3451 


1067 


6002 


1070 


6504 


1071 


5452 


1072 


2052 


1073 


3453 


1074 


6002 



110 



1075 

1076 

1077 

1100 

1101 

1102 

1103 

no4 

1105 

1106 

1107 

mo 

nn 

1112 

m3 

ni4 

ni5 

m6 

U17 

U20 



PROGRAM 5 (CONTINUED) 



Contents 

6514 

7500 

l4oi 

7600 

3460 

4024 

0101 

7000 

0031 

0022 

0024 

0030 

0101 

7000 

0041 

0030 

0002 

0030 

2030 

3060 



Comments 



Can A/D 



Read into A 

Subtract 'Old 
Value' 

Store E 1 

Multiply E 1 
by 100K 1 



Divide product 
by 100 



Load H 1 



Add 'Old Value' 
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PROGRAM 5 (CONTINUED) 



Location 


Contents 


Comments 


1121 


4030 


Store adjusted H 1 


1122 


7500 


Call D/A 


1123 


2401 




1124 


7303 


Output adjusted H 1 


1125 


0031 




1126 


7013 


Jump to second 
delay 


1127 


0030 


Positive Step 
Solution 


1150 


7500 


Call A/D 


1131 


1401 




1132 


7600 


Read into A 


1133 


3460 


Subtract ’Old 
Value' 


1134 


4024 


Store E° 


1135 


2424 


Complement E° 


1136 


4024 




1137 


0101 


Multiply complemented 
E° by 100 K° 



11 4o 


7000 


ll4l 


0031 


1142 


0020 


1143 


0024 


1144 


0030 



112 



1150 

1151 

1152 

1153 

1154 

1155 

1156 

1157 

11 6o 

ll6i 

1162 

1165 

11 64 

1165 

1166 

1167 

1170 



PROGRAM 5 (CONTINUED) 



Contents 


Comments 


0101 


Divide product 
"by 100 


7000 




004l 




0030 




0002 




0030 




2430 


Complement 

dividend 


4030 


Store H° 


2030 


Load H° 


3060 


Add 'Old Value' 


4030 


Store adjusted H° 


7500 


Call D/A 


2401 




7357 


Output adjusted H° 


0031 




2052 


Delay 


3452 




4052 




2050 




3450 
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PROGRAM 5 (CONTINUED) 



Location 


Contents 


Comments 


1171 


4050 




1172 


5450 




1173 


2050 




1174 


3451 




1175 


6002 




1176 


6504 




1177 


5452 




1200 


2502 




1201 


3453 




1202 


6002 




1203 


6514 




1204 


7500 


Call A/D 


1205 


l4oi 




1206 


7600 


Real into A 


1207 


346o 


Subtract *011 
Value* 


1210 


4024 


Store E 1 


1211 


2424 


Complement E 1 


1212 


4024 




1213 


0101 


Multiply complemented 
E 1 by 100K 1 


1214 


7000 




1215 


0031 
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PROGRAM 5 (CONTINUED) 



Location 


Contents 


Comments 


1216 


0022 




1217 


0024 





1220 


0030 


Divide product 
by 100 


1221 


0101 




1222 


7000 




1223 


0041 




1224 


0030 




1225 


0002 




1226 


0030 




1227 


2430 


Complement 

dividend 


1230 


4030 




1231 


2030 


Load H 1 


1232 


3060 


Add ’Old Value' 


1233 


4030 


Store adjusted H 1 


1234 


7500 


Call D/A 


1235 


2401 




1236 


7303 


Output adjusted H 1 


1237 


0031 




1240 


7013 


Jump to second 



delay 
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PROGRAM 5 (CONTINUED) 



Location 

1241 

1242 

1243 

1244 

1245 

1246 

1247 

1250 

1251 

1252 

1253 

1254 

1255 

1256 

1257 

1260 

1261 

1262 

1263 

1264 

1265 

1266 



Contents 

0030 

2052 

3452 
4052 
2050 

3450 
4050 

5450 

2050 

3451 
6002 
6504 
5452 
2052 

3453 

6002 

6514 

7500 

l4oi 

7600 

4o6o 

7010 



Comments 



Second Delay 



Bootstrap 
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PROGRAM 5 (CONTINUED) 



Location 


LOW CORE ALLOCATION 
Contents 


Comments 


0000 


.... 


Address of 


0002 


0144 


Subroutine ARITH 


0010 


1000 


Address of start 


0011 


1026 


of program 
Address of negative 


0012 


1130 


step solution 
Address of positive 


0013 


1164 


step solution 
Address of second 


0020-21 


— 


delay 

K° 


0022-23 


— 


K 1 


0024-25 


— 


Current error 


0030-31 


— 


Current H 


0051 


3403 


Preset total for 


0053 


0003 


0.1 second delay 
Number of passes 


0060 




through 0.1 second 
delay 

•Old Value* 


0062 


— 


'New Value' 


OO63 


— 


Increment 


OO65 


---- 


•Diff' 
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PROGRAM 6 



Due to the basic similarity of this program to the 
preceding programs, the comments presented here will be 
general in nature. 



Location 


Contents 


Comments 


1000 


7500 


SENSE LOOP 


1001 


l 4 oi 




1002 


7600 




1003 


4 o 64 




1004 


7500 




1005 


1402 




1006 


7600 




1007 


4065 




1010 


2064 


Sense error in 


1011 


3461 


System 1 


1012 


4067 




1013 


3063 




10l4 


6351 


If error, set 


1015 


2067 


E 1 Flag 


10 l 6 


3463 




1017 


6250 




1020 


2064 




1021 


4 o 6 i 




1022 


6101 




1023 


4o66 


If no error, D/A 


1024 


7500 


1025 


2401 


System 1 


1026 


7353 




1027 


0067 




1030 


2065 


Sense error in 


1031 


3462 


System 2 


1032 


4067 




1033 


3063 




1034 


6337 


If error, set 


1035 


2067 


E 2 Flag and 


1036 


3463 


jump to solution 


1037 


6236 


routine 


lo 4 o 


2065 




104i 


4o62 




1042 


6101 




1043 


4 o 66 




1044 


7500 
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1052 

1053 

1054 

1055 

1056 

1057 

lo6o 

1061 

1062 

1065 

1064 

1065 

1066 

1067 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1100 

1101 

1102 

1105 

1104 

1105 

1106 

1107 

1110 

1111 

1112 

1113 



PROGRAM 6 (CONTINUED) 



Contents 


Comments 


2402 


If no error. 


7333 


D/A System 2 


0067 


2 o 4 o 


If E 1 Flag set. 


6125 


jump to Solution 


6122 


Routine 


2050 




3450 


Basic Delay 


4050 




5450 




2050 




3451 




6002 




6504 




6463 


Return to start 


6564 


of Sense Loop 


5440 


6536 




20 4o 




3401 




4o4o 




6542 




5441 




7010 




2041 




3401 




4o4l 




7010 

0066 

0000 

0000 

0000 




2067 




3467 




4067 


SOLUTION ROUTINE 


2040 


E 1 Flag test 


6003 


6304 


If = 0, jump to 


6205 


E 2 Flag test 
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PROGRAM 6 (CONTINUED) 



Location 


Contents 


Comments 


1114 


7101 


If » -1, jump to 


1115 


1370 


Negative Step 


1116 


7101 


Solution for System 


1117 


1250 


Positive Step 


1120 


2042 


Solution for 


1121 


6 002 


System 1 


1122 


6137 


If J ^ 0, Jump to 


1123 


7500 


1124 


1401 




1125 


7600 


Compute and output 


1126 


3461 


H° (+) for System 1 


1127 


4020 




1130 


2420 




1131 


4020 




1132 


0101 




1135 


7000 




1134 


0031 




1135 


0020 




1136 


0022 




1157 


0016 




n4o 


0101 




ll4i 


7000 




11 42 


004l 




1143 


0016 




1144 


0004 




1145 


0016 




1146 


24l6 




1147 


3061 




1150 


4oi6 




1151 


7500 




1152 


2401 




U53 


7305 




1154 


0017 




1155 


5442 


J = J + 1 


1156 


7101 




1157 


1370 


Jump to E 2 Flag test 


1160 


0016 




1161 


5443 


1 = 1+1 


1162 


2043 
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PROGRAM 6 (CONTINUED) 



Location 



H65 

1164 

1165 

1166 

H67 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1200 

1201 

1202 

1203 

1204 

1205 

1206 
1207 

1210 

1211 

1212 

1213 

1214 

1215 

1216 
1217 



1220 

1221 

1222 

1223 

1224 

1225 
1226 
1227 

1230 

1231 



Contents 



3403 

6506 

2042 

3401 

6003 

7101 

1760 

7500 

l4oi 

7600 

3461 

4o66 

2466 

4o66 

3420 

6303 

7101 

1205 

0101 

7000 

0031 

0066 

0024 

0016 

0101 

7000 

0041 

0016 

0004 



0016 

2416 

3061 

4oi6 

7500 

2401 

7310 

0017 

5442 

2043 



Comments 



If I t (0003), 

Jump to E 2 Flag test 

If J f 1, jump to 
E X R Flag test 



Compute and output 
H 1 (+) for System 1 



J = J + 1 
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PROGRAM 6 (CONTINUED) 



Location 



1232 

1233 

1234 

1235 

1236 
1237 

1240 

1241 

1242 

1243 

1244 

1245 

1246 

1247 

1250 

1251 
1252 

1253 

1254 

1255 

1256 

1257 

1260 

1261 

1262 

1263 

1264 

1265 

1266 
1267 

1270 

1271 

1272 

1273 

1274 

1275 

1276 

1277 



Contents 


Comments 


3443 


1=0 


kokj , 




7101 


Junrp to E 2 Flag test 


1370 




0016 




0000 




0000 




0000 




0000 




0000 




0000 




0000 




0000 




0000 




2042 


Negative Step 


6002 


Solution for 


6135 


System 1 


7500 




l4oi 




7600 




3461 


(This solution is 


4020 


handled exactly as 




is its positive 


0101 


counterpart . Thus, 


7000 


no comments will he 


0031 


made.) 


0020 




0022 




0016 




0101 




7000 




0041 




0016 




0004 




0016 




2016 




3061 




4oi6 




7500 
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1301 

1302 

1303 

1304 

1305 

1306 

1307 

1310 

1311 

1312 

1313 

1314 

1315 

1316 

1317 

1320 

1321 

1322 

1323 

1324 

1325 

1326 

1327 

1330 

1331 

1332 

1333 

1334 

1335 

1336 

1337 

1340 

1341 

1342 

1343 

1344 

1345 

1346 

1347 



PROGRAM 6 (CONTINUED) 



Contents 



2401 

7305 

0017 

5442 
7101 
1370 
0016 

5443 

2043 

3403 

6506 

2042 

3401 

6003 

7101 

1760 

7500 

1401 

7600 

4o66 

3420 

6303 

7101 

1330 

0101 

7000 

0031 

0066 

0024 

0016 

0101 

7000 

Oo4i 

0016 

0004 

0016 

2016 

3061 

U016 

7500 



Comments 
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PROGRAM 6 (CONTINUED) 



Location 


Contents 


Comments 


1350 


2401 




1351 


7310 




1352 


0017 




1353 


5442 




135*1 


2043 




1355 


3443 




1356 


4043 




1357 


7101 




1360 


1370 




1361 


0016 




1362 


0000 




1363 


0000 




1364 


0000 




1365 


0000 




1366 


0000 




1367 


0000 




1370 


2041 


E 2 Flag test 


1371 


6003 




1372 


6304 


If = 0, jump to 


1373 


6205 


delay 


1374 


7101 




1375 


1650 


If = -1, jump to 


1376 


7101 


Negative Step 


1377 


1525 


Solution for 
System 2 


1400 


2044 




1401 


6002 


Positive Step 


1402 


6137 


Solution for 


1403 


7500 


System 2 


l 4 o 4 


1402 


If L 0. jump to 


1405 


7600 


1441 


l 4 o 6 


3462 




1407 


4030 


Compute and output 


1410 


2430 


H° (+) for System 2 


l4li 


4030 


1412 


0101 




1413 


7000 




l4l4 


0031 




1415 


0030 




l4i6 


0032 
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PROGRAM 6 (CONTINUED) 



Location 


Contents 


Comments 


1417 


0016 




1420 


0101 




1421 


7000 




1422 


004l 




1423 


0016 




1424 


0004 




1425 


0016 




1426 


2416 




1427 


3062 




1430 


4 oi 6 




1431 


7500 




1432 


2402 




1433 


7505 




1434 


0017 




1435 


5444 


L = L + 1 


1436 


7101 




1437 


1650 


Jump to delay 


l 44 o 


0016 




l44l 


5445 


K = K + 1 


1442 


2045 




1443 


3402 


If K ^ (0002), 


1444 


6506 


jump to delay 


1445 


2044 




1446 


3401 


If L ^ 1, jump to 


1447 


6003 


E 2 R Flag test 


1450 


7101 




1451 


2o4o 




1452 


7500 




1453 


1402 


Compute and output 


1454 


7600 


H 1 (+) for System 2 


1455 


3462 




1456 


4 o 66 




1457 


2466 




l 46 o 


3430 




1461 


6303 




1462 


7101 




1463 


1464 




1464 


0101 




1465 


7000 
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PROGRAM 6 (CONTINUED) 



Location 


Contents 


Conments 


1466 


0031 




1467 


0066 




1470 


0034 




1471 


0016 




1472 


0101 




1473 


7000 




1474 


0041 




1475 


0016 




1476 


0004 




1477 


0016 




1500 


24l6 




1501 


3062 




1502 


4016 




1503 


7500 




1504 


2402 




1505 


7310 




1506 


0017 




1507 


5444 


L = L + 1 


1510 


2045 




1511 


3445 


K = 0 


1512 


4o45 




1513 


7101 


Jump to delay 


1514 


1650 




1515 


0016 




1516 


0000 




1517 


0000 




1520 


0000 




1521 


0000 




1522 


0000 




1523 


0000 




1524 


0000 




1525 


2044 


Negative Step 


1526 


6002 


Solution for 


1527 


6135 


System 2 


1530 


7500 


(As before, no 


1531 


1402 


comments will 1 


1532 


7600 


for the negat: 


1533 


3462 


step solution. 


1534 


4030 
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1535 

1536 

1537 

1540 

1541 

1542 

1543 

1544 

1545 

1 546 

1547 

1550 

1551 

1552 

1553 

1554 

1555 

1556 

1557 

1560 

1561 

1562 

1563 

1564 

1565 

1566 

1567 

1570 

1571 

1572 

1573 

1574 

1575 

1576 

1577 

1600 

1601 

1602 

1603 



PROGRAM 6 (CONTINUED) 



Contents 



Comments 



0101 

7000 

0031 

0030 

0032 

0016 

0101 

7000 

004l 

0016 

0004 

0016 

2016 

3062 

4oi6 

7500 

2402 

7305 

0017 

5444 
7101 
1650 
0016 

5445 
2045 
3402 
6506 

2044 

3401 

6003 

7101 

2040 

7500 

1402 

7600 



3462 

4o66 

3430 

6303 
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1605 

1606 

1607 

1610 

1611 

1612 

1613 

l6i4 

1615 

1616 

1617 

1620 

1621 

1622 

1623 

162U 

1625 

1626 

1627 

1630 

1631 

1632 

1633 

1634 

1635 

1636 

1637 

l64o 

l64i 

1642 

1643 

1644 

1645 

1646 

1647 

1650 

1651 

1652 

1653 



PROGRAM 6 (CONTINUED) 



Contents 



Comments 



7101 

1606 

0101 

7000 

0031 

0066 

0034 

0016 

0101 

7000 

004l 

0016 

0004 

0016 

2016 

3062 

40l6 

7500 

2402 

7310 

0017 

5444 

2045 

3445 

4045 

7101 

1650 

0016 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

2050 Delay 

3450 

4050 

5450 
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PROGRAM 6 (CONTINUED) 



Location 


Contents 


Comments 


1654 


2050 




1655 


3451 




1656 


6002 




1657 


6504 




1660 


2040 


If E 1 Flag ^ 1 , 


1661 


6133 


sense error for 


1662 


7500 


System 1 


1663 


l4oi 




1664 


7600 




1665 


4o64 




1666 


3461 




1667 


4067 




1670 


3063 




1671 


6316 




1672 


2067 




1673 


3463 




1674 


6215 




1675 


2064 




1676 


4 o 6 i 




1677 


6101 




1700 


4o66 




1701 


7500 




1702 


2401 




1703 


7303 




1704 


0067 




1705 


6107 




1706 


0066 




1707 


5440 




1710 


6io4 




1711 


20 40 




1712 


3401 




1713 


4o4o 




1714 


2041 


If E 2 Flag = 0 , 


1715 


6133 


sense error for 


1716 


7500 


System 2 


1717 


1402 




1720 


7600 




1721 


4065 




1722 


3462 
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PROGRAM 6 (CONTINUED) 



Location 



1723 

1724 
1723 
1726 
1727 

1730 

1731 

1732 

1733 

1734 

1735 

1736 

1737 



1740 

1741 

1742 

1743 

1744 

1745 

1746 

1747 

1750 

1751 

1752 

1753 

1754 

1755 

1756 

1757 

1760 

1761 

1762 

1763 

1764 

1765 

1766 

1767 

1770 

1771 



Contents 



Comments 



4067 

3063 

6316 

2067 

3463 

6215 
2065 
4062 
6101 
40 66 

7500 

2402 

7303 



OO67 

7H0 

0066 

5441 

7105 

204l 

3401 

4o4l 



7010 

1105 

0000 

0000 

0000 

0000 

0000 

0000 



2046 

6132 

2042 

3442 

4042 

2043 

3443 

4043 



2040 

3440 



Jump to start of 
Solution Routine 



RETURN ROUTINE 
If E X R Flag £ 0, 
jump to ramp solution 
routine (not included) 

J = 0 

1=0 



E 1 Flag = 0 



1 30 



PROGRAM 6 (CONTINUED) 



Location 



1772 

1773 

1774 

1775 

1776 

1777 

2000 

2001 

2002 

2003 

2004 

2005 

2006 
2007 

2010 

2011 

2012 

2013 

2014 

2015 

2016 
2017 



2020 

2021 

2022 

2023 

2024 

2025 

2026 
2027 

2030 

2031 

2032 

2033 

2034 

2035 

2036 
2037 

2040 



Contents 



4o4o 

2041 

6112 

7500 

1401 

7600 

4061 

7500 

1402 
7600 

4062 
7011 
7500 

1401 

7600 

4o6l 

7010 

7101 

1105 

0000 

0000 

0000 



0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

2047 



Comments 



If E 2 Flag ^ 0 , 
bootstrap System 1 
and Jump to Solution 
Routine 



Bootstrap Systems 
1 and 2, and Jump 
to Sense Loop 
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PROGRAM 6 (CONTINUED) 



Location 


Contents 


Comments 


2o4l 


6132 


If E 2 R Flag f 0, 


2042 


2044 


jump to ramp 


2043 


3444 


solution routine 


2044 


4044 




2045 


2045 


L = 0 


2046 


3445 




2047 


4o45 


K = 0 


2050 


204i 


E 2 Flag = 0 


2051 


344i 




2052 


4o4i 




2053 


2040 


If E 1 Flag f 0, 


2054 


6112 


bootstrap System 2 


2055 


7500 


and jump to 


2056 


l4oi 


Solution Routine 


2057 


7600 




2060 


4o6i 




2061 


7500 


Bootstrap Systems 


2062 


1402 


1 and 2, and jump 


2063 


7600 


to Sense Loop 


2064 


4062 




2065 


7011 




2066 


7500 




2067 


1402 




2070 


7600 




2071 


4062 




2072 


7010 




2073 


7101 




2074 


4000 






LOW CORE ALLOCATION 




0000 




Address of Sub- 






routine ARITH 


0001 


0001 




0002 


— - 


Q 


0003 


.... 


P 


0004 


Ol44 




0010 


1105 




0011 


1000 





C 
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PROGRAM 6 (CONTINUED) 



Location 



Contents 



Comments 



0016-17 

0020-21 

0022-23 

0024-25 

0030-31 

0032-33 

0034-35 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 



H's for output 
E(l) 

100K°(1) 

100K X (1) 

E(2) 

100K°(2) 

100 ^( 2 ) 

E 1 Flag 

E 2 Flag 

J 

I 

L 

K 

E X R Flag 
E 2 R Flag 



0050 

0051 



0060 

006l 

0063 

0064 

0065 

0066 



0675 



0100 



Delay counter 
Preset total 
for 25 milli- 
second delay 

Old Value (l) 
Old Value (2) 
Increment 
New Value (l) 
New Value (2) 
Temp Storage 



Note: Cells 0070 through 0077 and cell 0007 are reserved for 

Subroutine ARITH. 
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